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ABSTRACT 


A software defined radio is a much more flexible platform than traditional, 
hardware implemented radios. By implementing radio functions in software, and putting 
those functions on a Field Programmable Gate Array (FPGA) chip, users will have the 
ability to download mission specific radio capabilities. This thesis examines a 
fundamental piece of the receiver, the Phase-Lock Loop (PLL), simulates a software 
PLL, and investigates the effects of fixed-point versus floating point mathematics 
required for an FPGA based PLL. With a fixed-point PLL simulator, figures of merit 
such as lock-time, lock range, and pull-in range are determined for typical signal-to-noise 
ratio (SNR) levels. 


V 



THIS PAGE INTENTIONALLY LEET BLANK 


VI 



TABLE OF CONTENTS 


L INTRODUCTION.1 

A. BACKGROUND.1 

1. The Phase-Lock Loop.2 

2, Fixed-Point Arithmetic.2 

B. OBJECTIVES.4 

C. RELATED RESEARCH.4 

1, Digital Delay Lock Loop.4 

2, Digital Downconversion and Channelization.5 

D. THESIS ORGANIZATION.5 

II, PHASE LOCK LOOP BASICS.7 

A. PRINCIPLES OF THE PHASE-LOCK LOOP.7 

1. Phase Detector (PD).8 

2. Loop Filter.9 


3, Voltage Controlled Oscillator / Numerically Controlled Oscillator.... 10 


B. OPERATION OF THE PHASE LOCK LOOP.11 

1, The Transfer Function.13 

a. Phase Step Applied to Input . 15 

b. Frequency Step Applied to Input. . 16 

c. Frequency Ramp Applied to Input . 1 7 

C. PERFORMANCE MEASURES OF THE PLL.17 

1. The Lock Range.18 

2. Lock Time.19 

3. The Pull-in Range.19 

4. Pull-in Time.19 

5. The Pull-out Range.19 

D. PERFORMANCE OF THE PLL IN NOISE.21 

III. FIXED-POINT ARITHMETIC.25 

A, FLOATING-POINT NUMBERS.25 

B. FIXED-POINT NUMBERS.26 

1, Fixed-point Number Representation.26 

2, Precision and Range in Fixed-point Arithmetic.28 

3, Errors in Fixed-point Numbers.29 

IV. FIXED-POINT PLL SIMULATION.31 

A, PROCEDURE FOR MODEL.31 

1, Floating-point MATLAB Simulation.31 

a. Input Signal Assumptions . 32 

b. Determination of PLL Parameters . 32 

c. The MATLAB Simulation . 35 

2, Fixed-point Simulink Model.38 

3, Fixed-point Model.39 

B, PERFORMANCE AND ANALYSIS OF SIMULINK PLL MODEL.41 











































1. Performance Measures of the Fixed-point Model.41 

a. Lock Range . 43 

b. Lock Time . 43 

c. Pull-in Range . 44 

d. Pull-in Time . 44 

e. Pull-out Range . 45 

f Frequency drift . 45 

2, Analysis of PLL.48 

a. Errors in the Output of the PLL . 48 

b. The 8-bit PLL . 54 

c. Effects of Changing PLL Parameters .55 

E, CONVERSION TO HARDWARE.55 

V. CONCLUSIONS.57 

A. CAPABILITIES AND LIMITATIONS.57 

B. RECOMMENDATIONS FOR FURTHER RESEARCH.58 

C. FINAL COMMENTS.59 

APPENDIX A.61 

LIST OF REFERENCES.69 

INITIAL DISTRIBUTION LIST.71 


viii 






















LIST OF FIGURES 


Figure 1. Ideal Software Defined Radio [3].xiii 

Figure 2. Simple Phase-Loek Loop.xv 

Figure 3. Phase-Loek Loop with Signals of Interest.7 

Figure 4. Bode Plot of Aetive PI Filter. (From Ref [1]).10 

Figure 5. Exeiting Funetions as Applied to Input of a PLL. (a) Phase Step.11 

(b) Frequeney Step, (e) Frequency Ramp. (From Ref [1]).11 

Figure 6. Block Diagram of PLL Transfer Function.14 

Figure 7. Time Response of PLL for different values of ^ and = 200.16 

Figure 8. Relationship of frequency ranges of a PLL. (From Ref [1].).20 

Figure 9. plotted as a function of SNR^, where is normalized to the natural 

frequency. (From Ref. [2].).22 

Figure 11. (a) Frequency Time Response of PLL. (b) Phase Error Time Response of 

PEE.37 

Eigure 12. Time Response of PEE for an Input Erequency Earger than the Eock Range. 

(a). Output of System Demonstrating the Signal Eocks on a Erequency of 

2350 Hz. (b). Plot of the Error of the System.38 

Eigure 13. (a). Plot of Eixed-point Model NCO Output, (b). Plot of Eixed-point Model 

Error.42 

Eigure 14. Vehicle Driving 60 mph Around a Tight Curve to Illustrate how the Doppler 

Shift can Cause a Erequency Ramp of a Signal.47 

Eigure 15. Error Plot of the Output of the Eixed-point NCO versus Time.53 

Eigure A1. Phase-Eock Eoop Simulink Model.64 

Eigure A2. Phase Detector Simulink Model.65 

Eigure A3. Simulink Active P-I Eoop Eilter.66 


IX 





















THIS PAGE INTENTIONALLY LEET BLANK 


X 



LIST OF TABLES 


Table 1. 


Range and Preeision of a 16-bit Fixed-point Data Type 


29 


XI 




THIS PAGE INTENTIONALLY LEET BLANK 



EXECUTIVE SUMMARY 


The explosion of wireless and PCS serviees over the last deeade has created 
numerous incompatible air interface standards. A subscriber to one service will find his 
phone rendered useless when roaming in the coverage area of another service, using a 
different standard. These competing transmission formats each might have their own 
unique modulation type, multiple access technique, error control methods, call set-up and 
handoff protocol and voice compression algorithms. The military has seen similar 
interoperability of radio standards between tactical radios used by coalition forces in 
Desert Storm. 

The need for a flexible communication platform, capable of interfacing with the 
numerous standards and formats has become apparent. A software radio, implementing 
traditional radio functions in software, gives the user this flexibility. The ideal software 
radio shown in Figure (1) would digitize the entire received signal spectrum using a high 
speed Analog-to-Digital Converter (ADC), perform all demodulation, data protocol and 
processing functions using a general-purpose digital signal processor (DSP). 



Figure 1. Ideal Software Defined Radio. (From Ref [3].) 
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Due to limitations on A/D conversion speed, sampling the direct RF spectrum is 
not typically an option. A practical software radio will usually incorporate an RF front 
end, which filters and downconverts a portion of the received signal spectrum to much 
lower IF frequency such as 21.4, 70 or 160 MHz. These frequencies can be digitized 
directly with current state of the art A/D converters. This technique is called “IF 
Sampling”. 

The digital signal processing (DSP) functions of the software radio can be 
implemented either with general-purpose DSP chips or with reconfigurable Field 
Programmable Gate Arrays (FPGAs). In either case, the processing functions are in the 
form of software, available for download to the DSP engine. The resulting unit would 
have the capability to be reconfigured for any radio signal format. This gives the user the 
ability to download mission specific radio requirements, using the same platform for 
numerous radio applications. 

This thesis will focus on the implementation issues for an important signal 
processing function common to most communications receivers; that is, the Phase-Lock 
Loop (PLL). The PLL can take several forms such as the Costas Loop for carrier 
recovery and tracking, the early-late gate synchronizer for baud timing recovery, and the 
delay-locked loop (DLL) for spreading sequence tracking in spread spectrum systems. In 
each case, the same PLL loop theory presented in this thesis applies. 

In short, the PLL is a feedback loop device, which locks onto a received signal, 
meaning it synchronizes its output in-phase and frequency with its input. The PLL can be 
broken down into its three component parts: 1.) the phase detector (PD), 2.) a loop 
lowpass filter, and 3.) a voltage controlled oscillator (VCO) or numerically controlled 
oscillator (NCO), the latter being used in a software PLL version. Figure 2.shows a 
simple PLL schematically. 
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Figure 2. Simple Phase-Lock Loop 

While a traditional analog receiver implements the PLL components exclusively 
in hardware, a software-defined radio requires these components to be available as 
downloadable software, using a digitized signal as the input. While a software simulation 
of a PLL has been a reality for some time, the advent of fast Field Programmable Gate 
Array (FPGA) technology makes this a useful concept for high data rate signals (those 
greater than roughly 2 Mbit/sec data rate). Currently, the software components can be 
programmed onto an FPGA and used real-time for DSP functions. 

When deciding on DSP implementations, one must consider whether fixed-point 
or floating point arithmetic and number representation will be used. The benefit of 
floating point is the large dynamic range associated with the floating-point number 
representations. When implementing algorithms in floating point, the designer typically 
does not have to worry about issues such as rounding or truncation error, or numeric 
overflow. The disadvantage of floating point is the increased computational resources 
required and the processing speed limitations. For this reason, floating point 
implementations of receiver algorithms on general purpose DSP are limited to relatively 
low data rate signals. 

One of the reasons FPGAs are fast enough to be used is the fact that they use 
fixed-point representations of numbers rather than floating point. Fixed-point 
representation is a much more efficient way for a computer to do arithmetic, because its 
essentially a binary representation of a decimal number. Hence, a computer can do 
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arithmetic at a much faster rate using fixed-point numbers. With the computationally 
heavy nature of any DSP application, the efficiencies of fixed-point arithmetic add up to 
significant savings in time. The trade off is that representing numbers in this way either 
reduces the range of the number, or its precision. As an example, using 8-bit fixed-point 
numbers, a range of-256 to 256 only has a precision of 2.0. For a precision of .01, the 
range of numbers available are -2 to 2. This constraint needs to be realized when 
programming using fixed-point arithmetic. 

The aim of this research is to model a software phase-lock loop, observe its 
performance, and convert it into a fixed-point implementation to determine the effects on 
important performance such as lock time and pull-in range. The questions that required 
answering were 

• How many bits are needed in a fixed-point implementation for acceptable 
performance? 

• What type of errors do a fixed-point implementation introduce into the 
output signal? 

• What effect does fixed-point arithmetic have on figures of merit of a PLL? 

• What kind of signal-to-noise ratios (SNR) are required to lock the PLL in 
an acceptable time 

The source code used initially was MATLAB, but to convert the model into fixed- 
point, the MATLAB extension Simulink was used. The final simulation model used 16- 
bit fixed-point arithmetic and locked with acceptable SNR’s, with only small errors to the 
output signal. Lock time and ranges were not changed, however a smaller pull-in range 
than the floating-point equivalent was encountered. Errors at the output of the PLL due 
to the fixed-point implementation are observed, but nothing of critical size. The fixed- 
point model performed comparable to an analog or floating-point model in all ways, and 
could used as a basis to build a PLL on a FPGA or DSP chip. 
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I. INTRODUCTION 


A, BACKGROUND 

The revolution in wireless teehnology in the last deeade has created a need for a fast, 
flexible, and light radio, compatible with the numerous transmission formats and standards. The 
concept of a software radio has been developed to fill that need. The ideal software radio 
digitizes the radio spectrum at the receive antenna, providing all demodulation, decryption, and 
signal processing in software. Such a tool would have the capability to download appropriate 
software, depending on the mission and the transmission format required. It could be 
reconfigured to accommodate any RF-band modulation or data format or transmission standard, 
capable of operating within any communications network. 

To realize the ideal software radio, the entire signal spectrum would be digitized at the 
antenna by an analog-to-digital converter (ADC). A practical SW radio will usually incorporate 
an RF front end which filters and downconverts a portion of the received signal spectrum to much 
lower IF frequency which can be digitized directly with current state of the art A/D converters. 
The digital signal processing (DSP) is then performed in software in a reprogrammable Field 
Programmable Gate Array (FPGA). Finally, the demodulated signal is sent through a digital-to- 
analog converter (DAC) to generate the audio or video output if required. The ideal software 
radio architecture is shown in Figure 1, where the DSP functions are done on an FPGA chip. The 
radio concept described above could be reprogrammed to accommodate any radio standard or air 
interface by downloading the appropriate software algorithm. 

The signal processing functions required by the software radio includes tuning, fdtering, 
demodulation and decryption. Tuning is accomplished by mixing the digitized signal with a 
digital local oscillator to down-convert the signal to baseband. This mixer is simply a signed 
multiplier, sample by sample in the digital implementation. Typically, a complex baseband signal 
representation is used. The mixer will output the in-phase and quadrature-phase components in 
sine and cosine waveforms. This mixing process creates unwanted frequencies, specifically at 
twice the receive frequency. This spectral component is removed using a finite impulse response 
filter (FIR). The digital filtering also bandlimits the samples to the bandwidth of the signal of 
interest (a process called “channelization”). 
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Demodulation is performed by a variety of methods, depending on the modulation type. 
This research assumes demodulation will be performed by phase locking to the incoming signal. 
An error signal is fed back to the tuner to adjust the local oscillator frequency in a closed loop 
system. This feedback system is simply a digital or software phase-lock loop (PLL). The 
implementation of such a PLL is the focus of this research. 

1, The Phase-Lock Loop 

A simple phase-lock loop is pictured in Figure 2 and can be seen to consist of three 
component parts: 1.) the phase detector, 2.) a loop fdter, and 3.) a voltage controlled oscillator 
(VCO) or numerically controlled oscillator (NCO). Its purpose is to lock on to the frequency and 
phase of the input signal. In this case, the purpose of the PLL is to create a phase coherent as a 
local oscillator in the receiver for demodulation. 

The phase detector or phase comparator compares the phase between the input signal and 
the output signal. It generates a signal proportional to the phase error, or difference between the 
two signals phase. This can be done numerous ways; the simplest is to multiply the two signals 
together. 

The loop filter block is a low-pass filter that removes the high frequency terms the come 
from the multiplication of the input and output signal, leaving only the phase error. 

The NCO takes the phase error from the loop filter output and adjusts its output sinusoidal 
signal to force the error to zero. This adjusted signal is the feedback signal that goes into the 
phase detector, producing a second phase error and the process repeats. 

This type of PLL is a second-order system. The error eventually settles to zero, but the 
output is a damped oscillation. This oscillation is governed by the parameters of the PLL. 


2. Fixed-Point Arithmetic 

Measurements of physical quantities can take on many numerical representations. For 
example, the number one thousand can be represented by 1000, 1E3, 10 , or one thousand. In this 

case, the same quantity is represented using four different syntaxes. Another example of the same 
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quantity, represented in different scales, is the boiling point of water. Water boils at 100 degrees 
Celsius, 212 degrees Fahrenheit, 373 degrees Kelvin, or 671.4 degrees Rankine. 

Determining an appropriate seale and representation depends on many faetors. Suppose 
you need to measure the speed of a vehiele. The numerieal values have a limited range for this 
applieation. The slowest it ean travel is 0 mph at a dead stop and the top speed of the vehiele has 
been determined to be 150 mph. If an 8-bit unsigned integer is required, values in the range of 0- 
255 are possible. 

A typieal approaeh would be to assign one bit per mile-per-hour, making the integer 0 a 
dead stop and 150 the vehiele’s top speed. This seheme, while easy to convert, neglects the use of 
the integers 151-255, wasting 40% of the number range. 

Another approach would be to set the integer zero to a dead stop and the integer 255 to the 
top speed. This scale gives us much greater precision, 0.58823 mph per bit, because all 256 
values are used in the vehicle speed. The conversion, however, requires a division of 1.7, a 
relatively expensive operation for fixed-point processors. 

The trend of recent technology is to implement control systems and digital signal 
processing functions on digital hardware. Within digital hardware, numbers are represented as 
either floating-point or fixed-point data types. The number of bits used to represent both data 
types is a fixed word size. The range of fixed-point representation is much smaller than for 
floating point, thus to avoid overflow and quantization errors, fixed-point representation must be 
scaled. If floating-point numbers can effectively represent real world values with virtually no 
error, why use a fixed-point based implementation? The answer, of course, is cost, size and 
processing speed. 

A fixed-point hardware platform is architecturally much simpler than its floating-point 
counterpart. This means cheaper manufacture of the product. In addition to manufacturing 
savings, if scaled properly, fixed-point arithmetic can be significantly faster, saving computation 
time. 

Because the logic circuits of fixed-point hardware are less complicated than that of 
floating point, the chip size can be much smaller, reducing power consumption. This means 
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smaller batteries, and reduced heat as a by-product, removing the need for an expensive and bulky 
heat sink. 

For a software radio, programmed on a FPGA where speed is a critical limitation, and size 
and battery life are a concern, a fixed-point implementation is the only plausible solution. 


B, OBJECTIVES 

This research focused on the effects of fixed-point arithmetic on a simulated PLL, 
specifically, what number of bits are required to successfully implement a PLL, how does this 
affect the figures of merit of a PLL such as pull-in range, lock time, and lock range, and to 
determine the signal-to-noise ratios (SNR) needed for a fixed-point PLL to be effective in a 
software radio. 

These objectives were accomplished using MATLAB’s Simulink software. A PLL 
simulation using only fixed-point arithmetic was implemented, taking as input a noisy signal. 
This simulation was done using 16-bit fixed-point arithmetic. Simulation results were consistent 
with theoretical results, with the fixed-point implementation simply adding “noise” to the signal, 
causing the output signal to have small errors due to the quantization effects of fixed-point 
arithmetic. 

C. RELATED RESEARCH 

The related research is in identifying the feasibility and requirements for other 
components to make a software radio realizable. Some of this research is being done concurrently 
with this research. The first is a digital delay lock loop (DDLL); the second is investigating 
digital downconversion and channelization. 

1. Digital Delay Lock Loop 

For a fully functional software radio, a real-time, spread-spectrum, signal-processing 
block is desired. An integral function of a direct sequence spread spectrum receiver is recovering 
the underlying narrowband data through a procedure called dispreading. The key component of 
dispreading is a DDLL. A DLL generates an exact replica of the spreading sequence generated at 
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the transmitter, allowing the reeovery of the narrowband data. This researeh was being eondueted 
by Captain Samuel Laboy, USMC. 

2. Digital Downconversion and Channelization 

A basie building bloek of any all-digital receiver is the downcoversion of a real signal at 
an intermediate frequency to complex in-phase and quadrature-phase components at baseband. 
The downcoversion component of the software radio consists of a complex NCO, local oscillators 
that act as mixers, and digital filter. Investigation in the feasibility of digital downconversion and 
performance characteristic compared to an analog equivalent is being pursued by Lieutenant 
Michael Snelling, USN. The aim of this research is to build a working MATLAB model and to 
analyze the results of the simulation. 


D, THESIS ORGANIZATION 

This thesis is organized to mirror the research. Chapter 2 is an overview of the operation 
of the PLL and derives the figures of merit for analysis of the final model. Chapter 3 is a tutorial 
on fixed-point arithmetic and explains the advantages that it offers over floating-point arithmetic. 
Chapter 4 is the culmination of this research. It begins by describing the building of a floating¬ 
point PLL and then the procedures for converting it to a fixed-point model in Simulink. Analysis 
of the performance of the fixed-point PLL according to the derived figures of merit derived in 
Chapter 2 is also done in Chapter 4. The conclusions for the feasibility of a FPGA based PLL and 
recommendations for future research are in Chapter 5. The Appendix gives the MATLAB code 
and Simulink model for the floating-point and fixed-point phase lock loops. 
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II. PHASE LOCK LOOP BASICS 


A, PRINCIPLES OF THE PHASE-LOCK LOOP 

A phase-lock loop is a circuit, or software, designed to track a given Reference signal in 
both frequency and phase. Its applications are far reaching from AM/FM radio demodulation, to 
television sets, to coders and decoders. The PLL has the ability to synchronize its output in-phase 
and frequency with an input signal, meaning the phase error between the PLL’s output and the 
input signal is zero, or remains constant. If a phase error is introduced, by a phase change or a 
frequency change, the PLL’s feedback control mechanism adjusts the oscillator’s output to 
account for it. 

To get a basic understanding of the PLL, the principles of the linear PLL will be 
examined, which is shown schematically in Figure 3. It consists of three blocks: the phase 
detector (PD), the loop filter (LF), and the VCO/NCO. For ease of Reference, the same notation 
and signal names used by [I] will be used here. Consequently the frequency of Uy{t) in 
radians/second is ft),, (O^ is the frequency of WjCOand 9^{t)is the phase error or difference in- 
phase of Wj(t) and Uj^{t) . With signals defined, a closer look at each functional block is in order. 



Figure 3. Phase-Lock Loop with Signals of Interest 
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1. Phase Detector (PD) 

The PD takes as its input two signals, the referenee signal, Wj(t), and the output of the 
VCO, ^2(0 • This block simply compares the phases of these two signals and produces a signal 
(t) proportional to the phase error 0^ , specifically, 

uAt) = KA ( 2 . 1 ) 

where is the PD gain in volts. Obtaining the signal Uj{t) can be as simple as a multiplier. 

The PD adopted in the simulation uses the in-phase and quadrature-phase portions of the input 
and output signal. To see how this works, lets assume a phase difference between the u^{t) and 
(t) of 0^. Defining the in-phase and quadrature-phase components as /; = t/j cos(<Wit), 

= f/j sin(<Wjt), /j = t /2 cos(<W 2 ^ + Qi - ^2 sin(<W 2 ^ + , so that 

Mj {t) = f/j (cos(<Wjt) -I- i sin(<Wjt)) 

W2 (0 = U2 (cos(<W2^ + + * sin(<W2t + )) 

where t/j and t/j are the amplitude of the respective signals and is related to the PD gain . 
Using the trigonometric identities 

cos Xcos y = ^[cos(.r -I- y) -I- cos(.r - y)] andsin vsin y = ^[cos(v - y) - cos(.r -I- y)] 

and a little algebra, 0^ can be extracted. For ease of computation, we assume 0)^ = 0 ) 2 . If this is 

not the case, the PD output has an additional frequency error, but this ac component will be 
filtered out by the loop filter explained below. 

I = I,^l2+Q,^Q2 

= t/j *U2[cos{cOit)cos{cjD2t + 0J + sm{cOit)sm{Q)2t + 0J] 

1 ( 2 . 2 ) 

= —U^*U2[cos(((O^+(O2)t + 0J + cos(0J + cos(0J-cos((a)^+(O2)t + 0J] 

= U,*U2COS(0J. 

Similarly, 
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Using equations (2.2) and (2.3), 


Q = h*Q2-Q,*h 


(2.3) 


0^ = tan ' 


I 1 ) 


As desired, the PD output signal (t) is the phase error. 


(2.4) 


2, Loop Filter 

As alluded to earlier, if the frequeneies of the two signals differ or alternate PD 
implementations are used, the output signal, u^{t), of the PD will have an unwanted ac 

component. This ac signal is superimposed on the desired dc component representing the phase 
error. To remove the ac component, a simple low pass loop filter is used. 

To implement a loop filter, numerous strategies exist. While a high-order finite impulse 
response (FIR) filter will ensure the ac component is removed, the delay in response is too great, 
where locking onto a signal as fast as possible is needed. A first-order low-pass filter has a quick 
response and successfully removes the unwanted oscillations. The most common loop filter used 
is called an active PI filter (PI = proportional + integral, taken from control theory and named due 
to the fact that it has a pole at s = 0, hence acts as an integrator) [1]. Taken from the 
corresponding RC circuit filters and implemented in software using their transfer functions, the 
filter has the Bode plot depicted in Figure 4 [1]. The transfer function for the PI filter is given by 

F(.) = l±^a (2.5) 

5Tj 

where in the analog circuit world, xi and X 2 are RC time constants of the circuit fdter and are 
determined by the loop bandwidth of the system, which is in turn a function of the noise. 
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Figure 4. Bode Plot of Active PI Filter. (From Ref [1]). 

3. Voltage Controlled Oscillator / Numerically Controlled Oscillator 

The NCO takes as its input the output from the loop filter, Uj- (t), which is proportional to 
the phase error. By using an appropriate scaling, or gain constant , of the NCO, we can adjust 
the output, MjCO? to account for the phase error in the two signals. The NCO has a center 
frequency. This is the oscillator frequency without any adjustments. This value is selected to 
reflect approximate working frequencies of the system. The PLL will properly lock onto a signal 
within a percentage range of the center frequency. Call this frequency roo- An adjustment made 
on the frequency of the NCO will be an offset from wq- Thus the output of the NCO has a 
frequency, 

ft>2 (0 = C0f^+ KqUj (t). (2.6) 

The discrete nature of the software system requires the output of the NCO to be a phase, 
as opposed to frequency, and thus by definition the phase is given by the integral over the 
frequency variation, 

= (2.7) 

Using a look-up table, the output of the system is taken as 

=cos(6>2(0) 

^li^QuadraturePhase ~ Sin(02 (0) 
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Either the in-phase or the quadrature-phase ean be used as the oseillator, but with the 
implementation deseribed above, both are needed as inputs to the PD. 


B. OPERATION OF THE PHASE LOCK LOOP 

To understand the operation of the PLL and how the bloeks work together, assume first 
that the system is initially loeked, i.e., , and examine three types of signals at the input; a 

phase step, a frequeney step, and a frequeney ramp. Taken from [1], the signals are shown in 
Figure 5. 





Figures. Exeiting Funetions as Applied to Input of a PFF. (a) Phase Step, 
(b) Frequeney Step, (e) Frequeney Ramp. (From Ref [1]). 


The aim is to eharaeterize these signals in terms of their phase, and proeeed with analysis 
using a generie phase term in the signal. We assume an input signal of a sinusoid as previously 
done, so that Mj(t) = f/j^ sin(&)jt-i-0j(t)), where 0^{t) is a generie phase signal that represents the 
various exeitation input signals. 

For the phase step, whieh represents a phase-modulated signal, 9y{t), simply performs a 
step ehange at t = 0, and is given by, 
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6>i(0 = A^>m(0 


( 2 . 9 ) 


where 0 < A^> < 2k is the magnitude of the phase ehange and u{t) is the step function. 

The next case is an example of a frequency-modulated signal. This signal has a step 
change in frequency at t = 0, which is given by the increment A<w. Thus the corresponding input 
signal is given by 

M, (t) = sin((<w, -I- AQ)tu{t))) 

= sin(<W]t -I- Acotu{t)). 

Thus 

e,{t) = Am. ( 2 . 10 ) 

This phase signal is just a ramp function for t > 0. 

The final signal is one whose frequency increases linearly with time. As shown in Figure 
5, the frequency is a ramp function, and the rate of increase in frequency is the slope of the line. 
If the slope of the ramp function for frequency is Ad), then the total frequency of Mj(t)is 
<W, -I- Am . By definition, the frequency of a signal is the first derivative of its phase with respect 
to time, 

0 ), +Ao) = — 
dt 

Thus the phase of the signal at time t is the integral of its angular frequency over the time interval 
0 < T < t, and so the input signal Mj (t) can be written as 


Mj (t) = t/j sin j(<yj -I- A(bT)dT 

0 


= t/j sin 


( 

(Od^-Aco 

V 



Consequently, the phase signal is given by 


ddt) = Aa— 


( 2 . 11 ) 
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1. The Transfer Functioni 


It is often helpful in examining the operation of any system to derive and understand the 
transfer funetion, H{s), of the system. This will be the approaeh in understanding the dynamie 
operation of the PLL. The transfer funetion is defined as the Laplaee transform of the output 
divided by the Laplaee transform of the input. For a PLL, we want to relate the phase signals of 
the output and input, thus 

= ^ ( 2 . 12 ) 

0,(s) 

where 0 , ( 5 ) is the Laplaee transform of the phase signal 6.{t) . To build this mathematieal model 
of the system, assume initially the PLL is loeked, so that <», = <^2 • This being the ease, the output 
of the PD is the phase error 6^ = 6^ times the PD gain eonstant , i.e.. 


whieh implies 


^d(t) = KA 


(2.13) 


(2.14)2 


The next bloek in the PLL is the loop filter. The transfer funetion of the PI aetive filter 
was already given in Equation (2.5). The last bloek to derive the transfer funetion is the 
VCO/NCO. Reeall that the NCO adjusted the output frequeney of u^it ),depending on its input 
from the loop fdter. A negative input redueed the output frequeney and a positive input inereased 
it. The eenter frequeney of the NCO, 0 )^, was the starting point for this adjustment. The output 

frequeney, (O^it ), was defined by 


<»2 (0 = <^0 + (0 


(2.15) 


where K^Uf (t) is the variation in the frequeney. However we want the output phase of the signal, 

not its angular frequeney. Using the definition of the phase as the integral over the frequeney 
variation, 

1 Transfer function analysis and derivation done using [1]. 

2 Note here we assume over a small time period, the phase error does not change dramatically and is thus 
approximated to be a constant. 
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( 2 . 16 ) 


^ 2(0 = Kq ^Uf{t)dt, 

The Laplace transform of an integral is 1/s, thus 

= (,) ^ = (2.17) 

s s 

The closed-loop block diagram of the transfer function is shown in Figure 6. 



Figure 6. Block Diagram of PLL Transfer Function 


From control theory, the closed loop transfer function for the entire system is given by 


_ K,K,F{s) 

\ + H,^{s)F{s)H^,^{s) s + K,K,F{sy 


(2.18) 


Simplifying this to standard notation. Equation (2.18) becomes 


H{s) = 


^ ° (1-|-5'T2) 


2 , 

5 ’ -\-s 


'kaF 


-I- 




(2.19) 


V 




J 


^2 


Control theorists put this in terms of natural frequency, and damping factor, ^, in 

order to use a standard equation which can easily be understood and studied. The equivalent 
transfer function equation in normalized form is given by 


H{s) = 


2sCo)„+0)l 

s^+2sCo)„+0)y 


( 2 . 20 ) 


14 





where 



^ 2 

To ealeulate the time response of the transfer function, the inverse Laplace transform of 
(2.20) is required. The time domain equivalent equation becomes 



which is simply a damped oscillation. 

To get a feel for the dynamic operation of the PLL for the three Reference inputs 
described above, we look at the transient response of the transfer function in normalized form for 
different values of (O^ and ^. Then finding a desired response, values to build the PLL can be 
determined. 

a. Phase Step Applied to Input 

To see a typical transient response, a phase step was applied to the input at t = 0, 
corresponding to the first type of input above, so that 

6>i(0 = A^>m(0 . 

Figure 7 shows the error time response of the transfer function versus time for various values of 
^ .3 A natural frequency of 200 rad/sec is used. The value of <W„only scales the x-axis, not the 

3 The time response can be obtained by taking the inverse Laplace transform of the transfer function. For this 
transfer function the time response is given by Equation (2.21) 
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response of the system or its oharaeteristies. To note is the fact that as the error 

approached zero. This can also be derived from the transfer function by finding the error transfer 
function, defined by 


liXs) = 


_ s 

0 ,( 5 ) s + K,K,F{s) 


( 2 . 22 ) 


and using the final value theorem of the Laplace transform which states 

0 ^( 00 ) = lim^^o 50^(5) = O. 



Figure 7. Time Response of PLL for different values of ^ and = 200 . 


b. Frequency Step Applied to Input. 

Again applying the final value theorem, a frequency step corresponds to a transfer 
function input of 


0 ,( 5 ) = Aco! s 
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Performing the same ealeulation on the phase error transfer function, 


0^(00) = lim^^o 


Ao) 

+ 2sCco„+co^„ 


Aco 


(o: 


(2.23) 


This value approaches zero for small frequency steps and high gain loops, which 
increases the value of 

c. Frequency Ramp Applied to Input 

Using similar techniques for a frequency ramp, the results are less encouraging. 

The error is 9^°°) = Ad)l co \, where Am is the rate of change of the input signal frequency. Thus 

for a large rate of change in frequency, the PLL unlocks. Experimen-tation with PLL’s has 
shown that a practical design limit to the rate of change in frequency is, 

=-^. (2.24) 


C. PERFORMANCE MEASURES OF THE PLL 

For a full derivation of the following key PLL parameters, one is encouraged to see 
Reference [1]. For the purposes of this research, it is sufficient to define the parameters and state 
the results for a given PLL, and shed some light on the competing parameters and their 
interaction. The focus will be on five key parameters that govern the dynamic performance of the 
PLL 

• The lock range, 

• The lock time, Tl 

• The pull-in range, A(Op 

• The pull-in time, Tp 

• The pull-out range AtOp^ . 

Additionally, a figure of merit for the PLL of loop bandwidth. Bp will be given 
consideration during the study of PLL performance in the presence of noise. 
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1. The Lock Range 

The lock range is defined as the range of frequency offset between the center frequency of 
the VCO and the reference frequency in which the PLL locks within one single-beat note between 
reference and output frequencies. We assume that the PLL is initially unlocked and switched on 
at time t = 0. To derive this figure of merit, we assume a simple PD of multiplying the input 
reference signal and the NCO output signal to produce the phase error. If this is the case, then for 
a frequency offset of A<y, 

(t) = sin(A<yt) -l- higher-frequency terms (2.25) 

where the higher-frequency terms are discarded due to the loop filter. Looking at the output of 
the loop fdter, the result, 

Uf (t) = I F{Aco) I sm(A<») (2.26) 

is obtained. The output of the LF, (t) , is just an ac signal causing a frequency modulation of 

the NCO, with a peak frequency deviation of \ F{A(o) \. For lock to occur, this peak 

frequency deviation must be less than the frequency offset Ary. If the offset is larger than this 
frequency deviation, lock cannot occur in a single cycle. The lock range can be determined by 
determining when this peak frequency deviation is just as large as the frequency offset. Thus the 
equation for locking becomes 

Aco,=K,K,\F{Aco,)\. (2.27) 

This is a non-linear equation, but to solve it a practical approximation can be made. From the 
Bode plot of the loop filter transfer function in Figure 4, the lock range is greater than the cut-off 
frequencies I/Tj and l/fj. Thus a conservative approximation is | F(Ary^) |= Tj / fj. Using 
Equation (2.20), the lock range becomes 

A(0^ - 2^ry„. (2.28) 

Simulations and experiments have shown that this is a conservative approximation and can be 
used confidently in the design process. Clearly a larger lock range is desirable. 
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1. Lock Time 


Lock time is the time the PLL takes to lock onto a signal when it is initially unlocked, 
given that the reference signal is within the lock range. When ^ < 1, the transients of the damped 
oscillation die out after one cycle, thus the lock time ean be approximated accurately as 



(2.29) 


2, The Pull-in Range 

The lock range is a subset of the pull-in range, whieh is defined as the range of frequeney 
offset in which the PLL will eventually lock after a number of cycles. This is dependant on the 
type of loop filter used. For the aetive PI filter, the pull-in range is infinite, i.e., 

>oo, (2.30) 

thus any referenee frequeney input to the PLL will eventually be locked onto. For different types 
of loop filters, with a finite gain at de, the pull-in range is deereased. See [1] for further details on 
pull-in range for different filter types. 


3. Pull-in Time 

Pull-in time is the time required for a PLL to look onto a frequeney within the pull-in 
range. Defining Aa)^ as the frequeney offset a)^-a) 2 , the pull-in time is determined to be 


AcDq 

16 C(^l 


(2.31) 


Note that, as Aco^ oo, the pull-in time approaches infinity. The larger the offset, the 
longer the pull-in time. 


4, The Pull-out Range 

The pull-out range is defined as that frequeney step which causes the PLL to unlock if 
applied to the reference input. Direot oaloulations of the pull-out range are not possible, but in 
[5], oomputer simulations have arrived at the following approximation 
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AO)po='^.SQ)SC + d- 


( 2 . 32 ) 


The relationship of the frequeney ranges is depieted in Figure 8. Starting from the eenter 
frequeney of the system, the loek-range is the smallest of the range parameters. Next is the pull¬ 
out range, and finally the pull-in range has the largest. The range of values outside the pull-in 
range is ealled dynamieally unstable due to the faet when operating dynamieally, the system will 
not look if a referenoe frequeney in that range is enoountered. Inside the pull-in range, and 
outside the look range is ealled oonditionally stable. This means the system will eventually look 
on frequenoies in this range, but oould unlook temporarily with a referenoe frequeney in this 
range. 



Figure 8. Relationship of frequeney ranges of a PLL. (From Ref. [1].) 
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D, PERFORMANCE OF THE PEL IN NOISE 


An in-depth analysis of the performanee of the PLL in noise is outside the seope of this 
researeh; however, some insight into what ean be expeeted from a PLL operating in the real world 
is required. The following definitions are used and assumptions made in this diseussion of noise. 

1 . SNRi=\0\og{Pg/P^), where P^ is the referenee signal power and is the noise 
power, is the signal-to-noise ratio at the input of the PLL. 

2. All noise is assumed to be Additive White Gaussian Noise (AWGN), meaning that it 
has a fiat power speetral density. 

3. A low-pass pre-filter of bandwidth B. is implemented prior to the PLL. 


4. The noise bandwidth is defined as = ^\H{jco)f Jtu, where//(jtu) is the transfer 

0 

funetion of the PLL. Inserting equation (2.20) for H{jco) , the integral solves to 


Bl= — 
2 


f 


r + — 


(2.33). 


5. The signal-to-noise ratio at the output is defined by 

= Sm, i. (2.34)4 

IB^ 

This figure is helpful in determining how often the PLL unloeks due to noise, and the 
SNR^ required to ensure a loek. 


Reeall that the loek range, A<w^, from Equation (2.28) is proportional to the natural 
frequeney and is desired to be as large as possible; henee, theoretieally, designing the PLL with a 
<w„as large as required is an easy task, but with the introduetion of noise, and Equation (2.33), the 
larger , the larger is 5^. This suggests a trade-off between redueing the noise suffieiently and 
a larger loek range, whieh is exaetly the ease. To inerease the loek range of the system, a larger 
SNR must be used. If only a small SNR is available, 5^ must be made small, whieh reduees 
A6>z.- 


4 For a derivation of SNR, one is direeted to Ref [1]. 
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The same is true of the loek time. This parameter is inversely proportional to co ^, and thus 
larger values of (O^ reduee . Again, by redueing is inereased. 

Praetical experiments with PLLs have shown that for stable operation, we need an 
SNRi^ > 6 dB [1], A system eould eventually loek with a lower SNRj ^, but the phase jitter would 
eause frequent unloeks and the system would be of little value. This leads to the final important 
figure, that of how often a PLL system, on average, will unloek. is defined as the average 

time interval between two unloeks of the system. Figure 9 is taken from [2], which depicts as 
a function of SNR^. For high SNR^, very difficult to find or measure, thus Figure 9 is 
provided to give an idea of the relationship. 



C 0.5 I'U 


SNH|^ 


Figure 9. plotted as a function of SNR ^, where is normalized to the natural frequency. 

(From Ref. [2].) 


A phase lock loop is a feedback system designed to lock onto the frequency and phase of 


the input reference signal. This is done using a phase detector, a loop filter, and a numerically 
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controlled oscillator. This chapter described the operation of these components and derived a 
transfer function for the system. Using this transfer function a typical time response for the 
system was plotted. Important figures of merit were defined and derived for the PLL. These 
performance measures are lock range, lock time, pull-in range, pull-in time and pull-out range. 
Finally a section describing the operation of the PLL in a noisy environment and how a noise can 
affect the figures of merit of a system was included. 

The objectives of this research are to implement the phase lock loop using fixed-point 
arithmetic. In order to do this a study of fixed-point arithmetic is needed. The following chapter 
introduces fixed-point arithmetic and some of the considerations that need to be looked at when 
using this type of representation. 
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III. FIXED-POINT ARITHMETIC 


Real world values are approximated to nearly any degree of error in eomputers using 
floating-point arithmetie. While the preeision afforded by floating point arithmetie is 
advantageous in numerous applieations, the speed at whieh a eomputer ean ealeulate and 
manipulate floating-point numbers is a disadvantage. Binary fixed-point arithmetie takes 
advantage of the proeessor shift instruetion when multiplying or dividing by two to speed up 
arithmetie, simply by representing numbers in a different fashion. The seeond advantage to fixed- 
point arithmetie is storage spaee. If the range of the real world values are known, and does not 
need 32-bits to eover the range adequately, fixed-point numbers ean be sealed to eover the range 
and use perhaps half number of bits. To gain an understanding of how fixed-point arithmetie 
works, an introduetion to floating-point representations is given. 


A, FLOATING-POINT NUMBERS 

The IEEE standard 754 has dominated most of today’s proeessors for floating point 
arithmetie. It speeifies four formats, the two most eommon being single-preeision and double- 
preeision. Single preeision uses 32-bits and double preeision uses 64-bits. Both double and 
single preeision formats eontain three eomponents; a sign bit (S), a fraetion field (f), and an 
exponent field (e). Eigure 10 shows how the bits are alloeated for IEEE standard 754 floating¬ 
point format for both single and double preeision. As the name implies, the fraetion field is fixed 
and the radix point, defined by the exponent field as deseribed below, is variable. Subsequently, 
for small integer number, most of the f field bits are zeros, henee wasted spaee. 

The exponent field is expressed as 2'’“'^^ for single preeision and for double 

preeision where e is a variable between 0 and 256. The field seales the fraetion field and plaees 
the radix point aeeordingly. The range of numbers able to be expressed 
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SINGLE PRECISIDV 




■'51 


(b) DDUBLE PREEISION 


Figure 10. IEEE 754 Eloating Point Number Eormat 

is a funetion of the exponent field, which ranges from 2“'^^ -10“^* (e = 0) to 2'^* ~10^* (e = 
256)5. 

The fraction field, as already alluded to, is fixed, and therefore is responsible for the 
precision of the number. The precision is defined as the distance between two subsequent 
representable numbers. This is simply 2“^ , where f is the number of bits in the f-field. Eor 
single precision floating point, the precision is 2”^^ ~ 10“’. 

B. FIXED-POINT NUMBERS 

As can probably be concluded, fixed-point numbers fix the radix point and the fraction 
field is the variable. This has two distinct advantages; 1.) by allowing the fraction to be variable 
in size, the storage space for an integer can be tailored to the range of values need, and 2.) 
arithmetic is essentially integer arithmetic and, when done in binary, this amounts to shifts, rather 
than cumbersome multiplications. 


1, Fixed-point Number Representation 

Fixed-point numbers can be specialized to unsigned and signed integers, or fractionals; 
however, the most general form will be considered here, which is a signed general fixed-point 
number, either an integer, fraction, or combination. To understand the usefulness of fixed-point 
numbers, a more intuitive look at decimal fixed-point arithmetic will be presented, followed by 

5 A bit is needed for the possibility of exceptional numbers such as infinity or NaN. 
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the extension to full binary fixed-point arithmetie. Fixed-point numbers in a proeessor are 
represented solely as integers. In order to properly insert the deeimal point, a proeessor has a 
syntax for all numbers. For instanee, if the proeessor uses a S.3.2 syntax, this means a sign bit 
(S), 3 integer digits, and 2 fraetional digits (realize that here a 16-bit range is being used to 
speeify the sign bit and the live deeimal digits). The number 238.15 in fixed-point S.3.2 notation 
is represented as the integer 23815. The original number is multiplied by 10 , where the exponent 
two is taken from the third field. To add two fixed-point numbers, eonversion to the notation of 
the number with the largest third field is required. A numerieal example will help to illuminate 
the proeess. Suppose we have the number above, 238.15, in S.3.2 notation and want to add the 
number 145.6 in S.4.1 notation. The proeessor has the integer numbers 23815 and 1456 stored, 
but to add the two, it must eonvert both numbers to S.3.2. The whole proeess is as follows 


238.15* 

23815 

100 


100 23815 10 1456 

-=- and 145.6* — =- 

100 100 10 10 

1456 23815 1456 10 38365 

10 “ 100 10 10 “ 100 


whieh is 383.65 in S.3.2 format. 


The payoff eomes when we extend this to binary fixed-point numbers. Instead of 
multiplying by powers of 10, we multiply by powers of 2. In a proeessor, this is aeeomplished 
with shifts. Shifts in a proeessor are extremely fast, thus speeding up eomputationally heavy 
applieations. 

Of eourse a proeessor does not keep traek of the format as deseribed above, but the aetions 
are equivalent. A fixed-point number is represented by 

V = SQ + B (3.1) 


where V is the real world number, B is the bias, whieh is used for signed numbers, Q is the 
integer that eneodes V, and S = 2^ is the slope or sealing of the number. The sealing determines 
where the radix point belongs. 
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2. Precision and Range in Fixed-point Arithmetic 

The range of a number gives the upper and lower limits of representation, while the 
preeision or resolution gives the distanee between two sueeessive numbers in the representation. 
In fixed-point numbers these are a funetion of the sealing, S, and the number of bits used or length 
of V. The range and preeision of a number are in eonstant eompetition. For a large range the 
resolution must be redueed. Conversely, for a high resolution the range is deereased. To 
understand this relationship, ealeulations of both of important parameters will be made. The 
terminology that is used will be a high or large resolution is desirable, whieh is equivalent to a 
large or high preeision. Often the term large preeision indieates that the distanee between two 
sueeessive representable numbers is large. For the purposes of this researeh, a large preeision 
indieates a high degree of aeeuraey, and a low or small preeision indieates a low degree of 
aeeuraey. 

Preeision is determined by the least signifieant bit (LSB), more speeifieally, what effeet 
ehanging the LSB from a zero to a one has on the number. Changing the LSB is the smallest 
inerement that ean be made, thus the effeet of ehanging the LSB determines the preeision. So 
what does the LSB represent? The sealing faetor, S, determines where the radix point belongs, 
henee, determines what the LSB represents. Suppose E = 0, i.e., 5 = 2°. This implies the radix 
point is to the right of the LSB, and thus a ehange from a zero to a one in the LSB results in a 
ehange in the number of 1. The resolution of this sealing is 1. For a better resolution, E is 
required to be negative. For E = -S, the resolution is 2~^ = 0.00390625. Generally, the 
preeision is defined as the sealing 5 = 2^. 

This inereased preeision, however, eomes at a priee. Assume we are using n bits to 
represent the fixed-point numbers. This means that a total of 2" possible numbers ean be 
represented. It should be elear that if a finite number of values ean be represented, the greater the 
preeision, the lesser the range. To determine the range, simply multiply the preeision by 2" 
whieh gives the total range of numbers available. Reeall the bias, B, is simply a sliding window 
for the range of numbers. Hardware represents numbers from 0 to 2", thus for signed numbers a 
negative bias, B = -2"“' is used, whieh slides the range between -2"“' and 2"“'. Table 1 shows 
the range and preeision of a 16-bit signed and unsigned fixed-point numbers for various sealings 

using a zero bias for the unsigned ease and B = 2"“' bias for the signed. 
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Scaling 

Precision 

Range of Signed 
Values (low,high) 

Range of Unsigned 
Values (low,high) 

2"^ 

0.125 

-4096,4095.875 

0,8191.875 

2-^ 

0.0625 

-2048,2047.9375 

0,4095.9375 

2"' 

0.03125 

-1024,1023.96875 

0,2047.96875 

2"' 

0.015625 

-512,511.984375 

0,1023.984375 

2"’ 

0.0078125 

-256,255.9921875 

0,511.9921875 

2"® 

0.00390625 

-128,127.99609375 

0,255.99609375 

2-10 

0.00097656 

-32,31.999023 

0,63.999023 

2-12 

0.000244140 

-8,7.9997559 

0,15.9997559 

2"'' 

0.000030518 

-1,0.99996948 

0,1.99996948 


Table 1. Range and Preeision of a 16-bit Fixed-point Data Type 


3. Errors in Fixed-point Numbers 

Due to the relationship between range and preeision in fixed-point arithmetie, errors oeeur 
between real world values and fixed-point representations of those values. Therefore, it is 
important to know the ranges and values of the required numbers and how important aeouraey is 
when using fixed-point arithmetie. Errors ean oeeur in essentially two ways, the value being 
represented is outside the range of the sealing used and error due to a too low resolution. 

The first kind of error is often easier to find. In an applieation, an out-of-range error, 
while not neeessarily raising an error message flag, ean be found through analysis. Most fixed- 
point proeessors, when eneountering a value outside the fixed-point range, saturate the fixed-point 
representation to the max or min of the range. This ean be a very bad approximation of the real 
world value. Thus knowing the range of values in an applieation is vital so that proper sealing 
ean be used. 
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A low resolution can also lead to errors which may be difficult to find. A processor 
representing a number between two successive fixed-point numbers does one of three things: it 
rounds down, rounds up, or rounds to the nearest representable value. To reduce errors as much 
as possible, rounding to the nearest representable value is preferred. This reduces the error to half 
the precision and was the approach used for the PLL simulation in Chapter IV. Care must be 
taken when arithmetic is done in fixed-point because adding, subtracting and multiplying can 
have dramatic effects on the size of numbers. Subtracting two similar numbers can result in a 
number smaller than the precision, hence equates to zero. Equally troubling is adding a small 
number to a large number. The larger number scaling must be used because the larger range is 
required, but this means low precision. If the small added number is smaller than the precision of 
the larger number, the addition will have no effect. It will be equivalent to adding zero. 

While increasing the number of bits can increase precision and range to any amount, this 
is not a valid solution for practicality reasons. In today’s commercial market, 32-bit fixed-point 
processors are the largest available. Application-specific DSPs are made with a larger number of 
bits, but these are not readily available nor cheap. Careful consideration of application values and 
ranges must therefore be a priority when using fixed-point mathematics. 


Using a fixed-point arithmetic implementation requires knowledge of the range of values 
and precision required for the application. The range and precision is determined by the scaling 
used. For the fixed-point PLL in the next chapter, each input and output range was determined 
from the floating-point model, then appropriate scaling used at each level. The next chapter 
details the procedure in developing the floating-point model and conversion to a fixed-point 
simulation. 
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IV. FIXED-POINT PEL SIMULATION 


The PLL simulation developed in this researeh used Mathwork’s Simulink software 
paekage. The reasons for ehoosing Simulink for this applieation are numerous. This software is 
designed for eontrol systems and has an imbedded fixed-point bloekset. This feature negated the 
requirement to generate eode to simulate fixed-point arithmetie. Another desirable feature of 
Simulink is the Real Time Workshop, whieh generates C eode from a simulation model. This 
feature will allow easy implementation on an FPGA as C ean be easily eonverted to VHDL, the 
language required for programming an FPGA. Finally, the Simulink software was available and 
familiar. There was no need to buy expensive modeling or simulation software, and beeause the 
software was familiar, learning a whole new language and set of proeedures was negated. 

A, PROCEDURE FOR MODEL 

The proeedure for building the model was first to eonstruet a floating-point simulation in 
MATLAB using the design from Referenee [1]. Flaving a working floating-point model gives 
expeeted results of a PLL simulation and a eomparison for the fixed-point model. The next step 
was to implement the floating-point model in Simulink, the bloek diagram based simulation 
software by Mathworks. This done, the Simulink model was eonverted to a fixed-point version 
where analysis ean begin. 

1. Floating-point MATLAB Simulation 

Any software simulation that takes a signal as input must rely on sampling that signal and 
using the diserete data points. Thus to the software, a veetor of data of a 20 hertz signal looks no 
different than a veetor of data bits of a 200 kHz signal. As long as the sampling frequeney is 
known and is larger than the Nyquist limit, the aetual frequeney of the input signal ean be 
arbitrary. The veetor representing the input signal is a veetor of numbers, regardless of the aetual 
frequeney of the analog signal it approximates. 
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a. Input signal Assumptions 

For the purposes of the simulation the center frequency, Fg, of the NCO was taken 
to be 2 kHz. This just means the expected frequency of the input signal is “around” 2 kHz. A 
sampling frequency , of 40 kHz was used. While a sampling frequency this high is not 

required for a 2 kHz signal, the resulting output is much easier to analyze visually. The input and 
output amplitudes of the signals were unity, with no DC offset. This assumption is justified by 
realizing that the input signal can simply be normalized by its amplitude if other than one, and 
any DC offset can be removed by filter prior to input into the PLL. A Hilbert transform is 
performed on the real input signal creating a complex signal with in-phase and quadrature-phase 
components. The signal was initially given a random frequency , within 2.5% of the center 
frequency of the NCO. This assured the input signal to be within the lock range. It was also 
given a random phase between 0 and 7t. The signal to noise ratio (SNR) at the input to the PLL 
was assumed to be SNR^ =15 dB. This noise level assured a lock was possible, but also gave an 
idea of how noise affected the system. A pre-filter , of 1 kHz around the expected input signal 
of 2 kHz was assumed. This would leave the input signal around 2 kHz undistorted, but reduce 
the effects of noise. Summarizing the input signal assumptions, 

• Fg = 2 kHz . 

• F^ = 20 kHz . 

• F,. = 1 kHz . 

• Wj (t) = cos{27rF^t -I- po) -i- i sin(2.;rFg,t -i- po) -i- n{t). 
where F^ g [1950,2050], po& [0,.7r], and n{t) is the noise with SNR^ =15 dB . 

b. Determination of PLL Parameters 

The first step in developing an algorithm for a software PLL is to determine 
appropriate parameters depending on the application and situation. Typical communications 
signals have an SNR of 15 dB, thus we are assured that the PLL will lock, although this will be 
verified once all parameters are determined. With the noise bandwidth, not a factor, 
determining the lock range is required. If the expected input signal has a frequency of 
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around 2 kHz, a fair estimate is that the actual frequency of the signal will not vary by more than 
5%. This means we need a lock range of 100 Hz, or > 200;r rad/s. Setting 
= 200;r rad/s and using equation (2.28), the natural frequency in terms of the damping factor 
is easily determined. 

rad/s. (4.1) 

From control theory, the transfer function of the system is optimally flat for f = 0.707, but from 
an analysis of noise, the noise bandwidth as a function of ^ is flattest for f = 0.5, thus a 
compromise of ^ = 0.6. With this information the resulting natural frequency is 
= 166.66.7r rad/s . 


Both PD gain and loop filter gain were set equal to one. The final parameters to be 
determined were the lead and lag constants of the loop filter. The transfer function of the loop 
filter was given in Equation (2.5). To determine the constants, the transfer function needs to be 
related to the time domain input-output. To do this, the transfer function of Equation (2.5) is 
converted to its equivalent z-transform. This is done using the bilinear transform, such that. 


where. 


F{z) = 


l + a,z”' 


a, = -1 


b 


1 + - 


tan 


r rj. \ 
^ 2^2 j 


(4.2)6 


6 For a complete discussion on bilinear transforms see [4]. 
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and T^=\l. 

The inverse of Equation (4.2) must be found in order to derive the eorresponding 
difference equation which will enable the loop filter to be implemented in software. Recall the 
transfer function in the z-domain, F{z), is the input divided by the output. From Figure 3, the 
input to the loop fdter in the sampled digital world is (n) and the output is Uj (n). Equation 
(4.2) can be written 

Uj.{z) = U,{z)Fiz). (4.3) 

Substituting in F{z), 


Uj. iz) = U, iz) ^ Uf (^)(1 + ) = U, {z){b, + b,z -^), (4.4) 

l + a,z 

and taking the inverse z-transform back to the time domain, the time difference equation becomes, 

Uf (n) = -a^Uj (n -1) + b^u^(n) + b^u^(n-\). (4.5) 

Equation (4.5) lends itself well to a software implementation. 

The topic of noise and whether a sufficient SNR at the input signal is present was 
glossed over previously. Because of the high noise tolerance of PEL’s, high SNR’s is more of an 
issue in decoding and demodulating the signal than it is in locking onto the signal. For this reason 
SNR’s for locking are exceeded to get better demodulation and bit error rate (BER) curves. 
However, this topic will be addressed now. 

Recall that the SNR at the input of the PEL, SNR^, is related to SNR^, the SNR at 

the output of the PEE by Equation (2.34). This relation is a function of two parameters: the 
prefilter bandwidth and the noise bandwidth. Assuming a typical prefilter bandwidth of 1000 Hz, 
the only variable in the relationship is the noise bandwidth, . The noise bandwidth is related to 

the natural frequency and the damping factor by Equation (2.33). This equation results in 
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= 63.5 Hz . Equation (2.34) gives the output SNR^ = 23.9 dB. Experimentation has shown 
that fee’s will look with SNR^ > 6 dB, thus the PEE with the above parameters should be well 
with the looking limits with regards to noise. The simulation parameters found are listed below 
for ease of referenoe. 

• =166.66.7r rad/s. 

• f = 0.6. 

• Ao)^ = IOOtt rad/s (100 Hz) . 

• T^= 3.64e~\ = 2.3e~^. 

• a, =-l, ho =635.1475, h; =-621.4397. 


c. The MATLAB Simulation 

With the parameters established above, the algorithm needed to be determined. 
The feedbaok loop of the PEL was done using a simple MATLAB for-end loop. This was done 
for oonvenienoe beoause an input signal of a given length was assumed for the simulation, thus 
the number of times required to “feedbaok” was known. This simulated the faot that real radio 
signal will be sent in “blooks” or “paokets” in whioh the length of eaoh paoket is known. In s ide 
the loop are the three blooks of the PEL. 

The first blook is the phase deteotor (PD). This blook was implemented exaotly the 
way the PD in Chapter 2 was desoribed by using the in-phase and quadrature-phase oomponents 
of the signal and some trigonometrio properties. Instead of oontinuous signals being multiplied 
and manipulated, the sampled data at oorresponding time intervals was manipulated. The output, 
(n) , of the PD is the soalar number representing the phase error at time n. 

The signal (n) was then sent to the loop filter to get rid of the higher frequenoy 
terms from the PD. The loop filter was implemented direotly from Equation (4.5), produoing the 
output signal Uj- (n ), whioh is the error in the two signals at the PD. 

The final blook in the loop is the NCO. This blook needed to take the output signal 
of the loop filter and make a oorreotion to its osoillation depending on the magnitude of (n). 
From the disoussion of the NCO in Chapter 2.A.3, the phase of the oontinuous time signal is 
given by Equation (2.7). In the disorete realm, the phase change of ^ 2(0 can be determined by 
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( 4 . 6 ) 


A 02 — (<^0 (^))^5 • 

where 0)^ is the eenter frequeney of the NCO in rad/s. With the ehange in phase known, the total 
phase of the signal ean be determined by 

6>2(n +1) = 6»2 (n) + ((Do + K^u^ (n))!^ . (4.7). 


If 02 is initializd to 02(O) = O 5 this computation beeomes possible. Notiee that Equation (4.7) is 
iterative and theoretieally 02 can continue to grow, thus a simple if-then statement is used to 
bound 02 between -tt and tt . This is done by subtracting Itt from 02 whenever 02>7t. 
Beeause the output of the NCO is required to be the in-phase and quadrature-phase eomponents, a 
look-up table is used to get the desired output. 


U2(n)2 = cos(02(n)) 
U2(n)Q = sin(02(^))- 


(4.8). 


This signal is now sent to the PD, and the loop is eomplete. A eomplete listing of the eode used 
for the MATLAB simulation is given in Appendix A. 

While a full analysis of the performanee of the above PEL is not the goal nor 
helpful in this researeh, eonfirming the PEE algorithm operated eorreetly is required. The input 
signal used had the assumptions deseribed in Subseetion 4.A.l.a. However, instead of a random 
frequeney, an input signal of 2100 Hz is used. This signal represents a 5% error in the eenter 
frequeney of the model and the input signal, whieh is the worst ease expeeted. Eigure 11 shows 
the results of the simulation. The output frequeney of the NCO is shown in (a) and the phase 
error of the two signals is shown in (b). The output ean be seen to oscillate, but is damped and 
settles on the eorreet frequeney. The error plot ean also be seen to oseillate, but dies out and goes 
to zero. 

Evaluating look-time, reoall, = 2.;r/<y„ = 0.0120 s, whieh from the plots in 

Eigure 11 appear to be aoourate. The oscillations have died out at this time and look has been 
aohieved. 
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Figure 11. (a) Frequency Time Response of PLL. (b) Phase Error Time Response of PEL. 

Similar evaluations on the pull-in range, which is theoretically infinite, proved to 
be as accurate as testable in the simulation. The pull-in time of the system is exponential with the 
initial frequency displacement. Proving this is the case is not plausible for all ranges of 
frequencies. Flowever, given an initial frequency larger than the lock range of the system, the 
PEE indeed pulled the NCO output to the larger frequency offset and eventually did achieve lock. 
As seen in Eigure 12, the output of the PD eventually reaches zero after going through several 
oscillations. This figure was generated with an input signal of 2350 Flz, which is an offset from 
the center frequency of the system of 350 Flz, clearly outside the lock range of 100 Hz. From 
Equation (2.31), the pull-in time of the system should be 0.0346 seconds. This is a conservative 
time estimate. This can be verified from the plot, where the oscillations are completely damped 
out. 
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Figure 12. Time Response of PLL for an Input Frequency Larger than the Lock Range, (a). 

Output of System Demonstrating the Signal Locks on a Frequency of 2350 FIz. (b). Plot of the 

Error of the System. 

This concluded the floating point MATLAB simulation. The next step was to 
build this system in the control simulation language of Simulink. 

2, Fixed-Point Simulink Model 

The fixed-point blockset in MATLAB is a toolbox used in Simuli nk . Simulink is a 
supplement to MATLAB for building models and running simulations, thus ideal for the purposes 
of the PLL. It uses blocks as functions and mathematical operators. Implementing the functional 
blocks of the PLL was a matter of translating the MATLAB code to blocks in Simulink. The 
model can be seen in Appendix A. 

The Simulink model was built using the same parameters as the MATLAB simulation. As 
expected, its performance was identical to the performance of the MATLAB coded model, thus 
no results will be provided here. 
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3. Fixed-Point Model 


The point of this thesis was to establish the plausibility of a fixed-point PLL in software 
and to determine what type of performance can be expected from such a system. This was to be 
achieved by building a model and running simulations to analyze the results. In building the 
fixed-point model, the same PLL parameters as the floating-point model were used. 

As will be the case throughout the fixed-point process, the design problem faced at the 
input and output of each block is determining an appropriate scaling for the fixed-point number so 
that a large enough range is used, without introducing unacceptable inaccuracies due to a lower 
precision. As will be explained in Chapter 4.C, a 16-bit fixed-point implementation was used. 
Essentially, the reason for this was while an 11-bit representation was sufficient to run correctly, 
typical FPGA’s and DSP chips use 8-bit, 16-bit, or 32-bit fixed-point arithmetic. The 8-bit 
representation didn’t give enough precision for the system to lock with any reliability, thus 16-bits 
were used. 

The procedure to convert the floating-point Simulink model to a fixed-point model 
involved looking at each block and determining the range of the signal at the input of the block 
and the range of the signal at the output of the block. With this information, a scaling appropriate 
for the range can be determined and, subsequently, the precision can be found at that scaling. If 
the precision was too large, an alternate method must be used, or an increase in the number of bits 
used. 

To begin, the input signal to the PLL needed to be converted to its fixed-point equivalent. 
Recall the input signal was normalized to ±1. From Table 1, the best fit for this range of numbers 
is using a scaling of 2“'^. Because the range at this scaling is exactly the expected range of the 
input, any noise could be amplified, hence a slightly larger range was used, specifically a scaling 
of 2“'"^. This gives a precision of 6.1035e“% which is still below the noise of the input signal, 
thus is acceptable resolution. 

With the input signal converted to fixed-point, the logical next step was to convert the PD 

to a fixed-point equivalent. The block diagram of the PD is shown in the Appendix in Figure A2 

(duplicated on the next page for convenience). Tracing the connections, one can see that this is 

equivalent to the MATLAB code implementation in Appendix A. The PD takes the input signal 
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and multiplies it by the in-phase and quadrature-phase components of the NCO. In the block 
diagram of Figure A2, In2 (short for input 2) is the phase, 02 W > which goes into a look-up table 
for sine and cosine producing the in-phase and quadrature-phase signals. The multiplication 
mentioned above multiplies two signals with a range of values between -1 and +1, thus the range 
of the output of this block must be within the range of-1 and +1. As can be seen in Figure A2, 
below the multiplication blocks Productl through Productd, one can see the scaling is again 2“'"^, 
reflecting the range required at the output of this block. 
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Figure A2. Phase Detector Simulink Model. (Duplicated From the Appendix for Ease of 

Reference.) 

To belabor the point, the next blocks in the PD, Suml and Sum2, simply add or subtract the 
output of the multiplication blocks, producing a range of -2 to +2, which is the range of the 
scaling 2“'"^. As can be surmised, this process of determining the range at every block needs to be 
done, and appropriate scaling of the output determined. Care must be taken in ensuring that when 
applying the best possible range to the output, a resolution problem isn’t introduced. This issue 
will be discussed in Section 4.C. and is the reason an 8-bit implementation was not possible. The 
full block diagram of the Simulink PEL is included in Appendix A, with scaling at each block 
indicated. 
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The LF was designed using the 1/z block. Referring back to the z-transform of the 
difference equation established for the MATLAB simulation and rearranging equation (4.4), 

Uf (z) = -a,z~'U ^ (z) + hot/^z) + b,z~'U ^ (z), (4.9) 

where z~' is the unit delay operator. This equation is the one used in the Simulink model for the 
loop filter. 

The NCO was implemented similarly using constant blocks and multiplier blocks, 
ensuring the range of the inputs and outputs were scaled properly. The complicated circuitry in 
the NCO is due to the fact that 0^ grows each increment and must be reduced to -7t to 7t . This 
was done using a comparator block and a switch. The output of the switch is the top input, which 
is simply 9^ for ^2 < • If ^2 > ’ the comparator resulted in zero, the switch outputs its third 

input vice its first, which is 9^-ln . 

B, PERFORMANCE AND ANALYSIS OF SIMULINK PLL MODEL 

The performance of the model was consistent with the MATLAB simulation considered in 
Subsection 4.A.I.C. For high SNR values, the fixed-point model performed identically to the 
floating-point model in lock range, lock time and pull-in time. The error of the output signal due 
to the quantization of the fixed-point signal became evident in the fixed-point model. This error in 
effect reduces the SNR of the output signal. 

1, Performance Measures of the Fixed-point Model 

Figure 13 shows the time response of the output of the PD, corresponding to the error in 
the PLL. It is clear that the error, 6 ^, dies to zero, thus confirming the PLL achieves lock. To 
generate this plot a reference frequency of 2100 Hz was used. 
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Figure 13. (a). Plot of Fixed-point Model NCO Output, (b). Plot of Fixed-point Model Error. 

To establish the performance of the PEL, the theoretical values of the performance 
measures will be computed and then compared to the results obtained by the simulation. Erom 
Section 2.C, five key performance measures were derived. These are lock-range, lock-time, pull- 
in range, pull-in time, and pull-out range. Eor the parameters used in the simulation, the 
following theoretical values for the fixed-point PEL are calculated. 

1. = 200.;r rad/s. 

2. =0.0120 s. 

3. AcOp = oo. 


4. 


AO)q 


5. AcOpQ =1508 rad/s. 


Eor the testing of the fixed-point PEL, a worst-case scenario input was adopted. With this 
in mind, the random phase was dropped in favor of an artificially placed phase offset of 90 
degrees. This represents the worst possible phase offset and, thus, the most conservative 
measures of performance. 


42 








a. Lock Range 

The lock range of the PLL model is the range of frequencies at which the system 
locks in a single operating cycle. This corresponds to one period of the natural frequency . To 

determine this value experimentally, the input signal frequency was iteratively increased until the 
system did not lock in a single cycle. This established a bracket, or upper and lower bound on the 
lock range. By successively narrowing this bracket down, the lock range was determined. Using 
this procedure the lock range was found, Ao)^ =5\%7r rad/s, a value over twice as large as the 
theoretical value of 200.^ rad/s. 

While a large lock range is desirable, some explanation as to why it is larger than 
predicted is required. There are two reasons this is the case. First, when determining the lock 
range, an approximation on the gain of the loop filter at the frequency offset was used. The 
assumption was that the lock range was larger than the comer frequencies of the two time 
constants, and thus the comer frequency was used as an approximation. The lock range is indeed 
larger than the corner frequencies, but with the non-linearity of the loop filter transfer function, a 
tighter bound is difficult to calculate, thus the overly conservative loop filter comer frequency is 
used. 

Second, the fact that noise was not a factor in the design of this PLL, the noise 
bandwidth is very small, which means the noise is not hampering the locking of the device. This 
results in a larger lock range than theoretically calculated. 

b. Lock Time 

The lock time can be visually seen in the plot of Figure 13, but for a more accurate 
determination, the discrete output values were analyzed and the lock-time determined to be when 
the output entered and remained within 2% of the final value of the system. This is an accepted 
percentage for the settling time of a system in control theory. This is equivalent to the time 
required for the error to die down and remain below 0.02. However, the actual analysis is more 
involved. With the introduction of noise, a large noise value at any given sample can throw the 
output outside the 2% range, even though the output due to the system is within the 2% range. 
Examining the numerical values at each sample and finding where the output stayed below 0.02, 
the lock-time was found to be at 0.0128 s. Comparing this value to the theoretical value of 0.0120 
s, the fixed-point representation does not degrade performance in this area significantly. It must 
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be understood also that the theoretieal value of 0.0120 is assuming no noise, and that any 
signifieant spike in the noise at the output of the system eould effeet the settling time of the 
system when analyzing the output at eaeh time sample. However, it should be noted that running 
the simulation several times, the above settling time of 0.0128 was the average of all simulations. 

c. Pull-in Range 

The pull-in time for the seeond-order PLL is theoretically infinite, but an analysis 
of the pull-in range of the fixed-point system needs to be done. At some point, the frequency 
offset will get so large that the fixed-point precision or range will be exceeded and critical errors 
will be introduced. Incrementally increasing the input frequency, the fixed-point PLL eventually 
obtained lock for a max frequency of 2830 Hz. This value represents a 41.5% error from the 
center frequency of the PLL. 

The limitation seems to come from the loop filter. The range or precision, 
depending on the scaling used, is exceeded at the output of the loop filter, causing the input to the 
NCO to be inaccurate. The same limitation was discovered when attempting to develop an 8-bit 
fixed-point PLL and will be discussed in more detail in Section 4.D. 

d. Pull-in Time 

The pull-in time, Tp , of the system is dependant on the frequency offset from the 
center frequency. Realize that if the frequency offset is inside the lock range, the lock time is 
used and is the same for any frequency offset, as long as it is smaller than the lock range. The 
pull-in time varies with the magnitude of the offset. To get an idea of the performance as related 
to pull-in time, five input frequencies were used and compared to the theoretical values for those 
frequency offsets. The settling time again was calculated using a 2% tolerance. Table 2 shows 
the results for the five frequencies. The pull-in time is given in seconds for the theoretical as well 
as the experimental and the frequency offset is given in Hz. An initial offset of 300 Hz was 
chosen due to the fact that from the analysis of lock-range, = 259 Hz and, thus, a larger 
value than this was needed. The pull-in range of the system was determined to be 830 Hz, thus 
only values up to 800 Hz were included. 


Aw [Hz] 

Tp [s] (theoretical) 

Tp [s] 



(experimental) 
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0.0254 



0.0452 


500 

0.0707 

0.0446 

600 

0.1018 


700 

0.1385 

0.0836 

800 

0.1810 

0.1114 


Table 2. Pull-in Times for Various Frequeney Offsets 


The disagreement between the theoretieal and experimental values, espeeially as 
the frequeney offset inereases, eannot entirely be explained by the approximations used in 
deriving Tp . Another explanation lies in the diserete quantization of the fixed-point PLL. The 
output of the PD, or the error of the system, has a certain precision. The range of numbers 
required at that block, as described previously, determines this precision. If the value at the 
output of the PD is smaller than this precision, the value is rounded to zero. By the nature of the 
second-order PLL, any time two successive outputs fall within the precision limits, the output 
remains at zero. The above theoretical values assume an analog PLL, where precision is infinite. 
A fixed-point software PLL has quantization error, which in this case, output a zero even though 
the real-world value is larger than that. 

e. Pull-out Range 

The pull-out range is an important figure because it allows the user to determine 
how large a frequency step can be supplied to the input without unlocking the system. This 
assumes the system is already in lock and a frequency step occurs, as in an FM system. This 
value was tested by giving the PLL a signal, allowing lock to occur, then stepping up the 
frequency until it unlocked. Simulations of the fixed-point PLL agreed with the expected figure, 
which was obtained through computer modelling in Reference [5]. The fixed-point PLL 
experimentally had a pull-out range of 1501 rad/s or 239 Hz, compared to the theoretical value of 
1508 rad/s, or 240 Hz. Recall that the lock-range of the system was 259 Hz, thus if a frequency 
step larger than the pull-out range, but less than the lock range, the system will again lock in the 
lock-time, or 0.0120 seconds for this system. If the frequency step is larger than the lock-range, 
the system will acquire lock using the longer pull-in time as described above. 

/. Frequency drift 
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The magnitude of a frequeney drift is of interest beeause, for a software radio 
implementation and downloadable software radio packages, satellite communication is required. 
As neither satellites nor users are stationary, a frequency drift can occur due to the Doppler shift 
in frequency. The Doppler shift is governed by two parameters: the velocity and direction of 
travel relative to the antenna, and the wavelength of the signal. Specifically, the Doppler shift in 
frequency is, 

f,=jcose. (4.10) 

where v = velocity, A = wavelength or d f^, c=speed of light, and 0 = direction of travel from 
the antenna. For a frequency drift to occur, the satellite would need to either be changing velocity 
or changing direction. With an orbital satellite, the direction of movement is constantly changing. 
Therefore, an idea of what kind of drift the PLL can handle is an interesting figure of merit. 

From Subsection 2.B.1.C, the rate of change in a frequency ramp that will cause a 
system to unlock is governed by Equation (2.24). To test the fixed-point PLL model, a frequency 
ramp was applied to the input after the system was allowed to lock. The slope of this ramp was 
increased in each run until the system was observed to unlock and to be unable to track the 
frequency change. The rate at which the PLL was unable to track the frequency ramp was much 
smaller than that predicted by Equation (2.24), but this value was assuming no noise and an 
analog system. Because of the non-linearity of a frequency ramp, noise and quantization error are 
magnified. The results of the fixed-point PLL are =15,000, compared to the theoretical 

value of 130,000. The PLL could maintain lock up to this rate of change for about 0.25 seconds, 
corresponding to about 30 cycles of the system. Lrom 0.25 seconds on, the large frequencies 
were out of range of the fixed-point representation. The PLL maintains lock longer for a more 
gradual frequency ramp slope (see example below). The limiting factor is simply the range of 
frequency that must be covered for the case of a frequency drift. Lor a fixed-point PLL with a 
center frequency of 2000 Hz, the system unlocked when the reference frequency reached around 
3300 Hz. Lor smaller sloped frequency ramps, the system remained locked longer, but once it 
reached the upper bound of 3300 Hz, the system fell apart. 

To understand if this is a limitation for a real world application, an example of a 
moving user is given. There are two ways a frequency ramp can be seen at the receiver of a PLL 
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due to the movement of a user. The first is if the user is accelerating in speed. The second is if 
the user is changing direction. The second case will be considered first. The largest change in 
frequency due to a change in direction occurs when 0 goes from 1 to -1. This represents a user 
turning around and moving in the opposite direction. The user will be assumed to travel at a 
constant velocity of 60 mph (26.82 m/s). If the user is assumed to turn-around on a dime, the 
result is a frequency step. The performance of a frequency step was covered in Subsection 
d.B.l.e, therefore a circular path of the user will be assumed, which represents a constant change 
in angle, thus a frequency ramp as the input signal of the PLL. As shown in Figure 14, the user 
makes a tight circle with a radius of 10 m. To get anything but negligible values, a high 
frequency signal of 100 MHz is used. At 26.82 m/s the turnaround would take them 2.34 
seconds. The value 0 in Equation (4.10) goes from 0 to .tt in 2.34 seconds. This represents a 
slope of the frequency ramp of 48 radians. A slope of 48 radians would take 170 seconds, or 
almost three minutes for the PLL to unlock. The frequency ramp is only applied for 2.34 seconds. 
The system can easily handle a frequency ramp of this slope, even with an elevated frequency of 
100 MHz. 



Ligure 14. Vehicle Driving 60 mph Around a Tight Curve to Illustrate how the Doppler Shift can 

Cause a Lrequency Ramp of a Signal. 

The second case is when the user is accelerating. Suppose the user goes from a 
dead stop to 60 mph in 3.1 seconds (the fastest production car ever built is the Porsche Carrera 
twin turbo which was tested at 0-60 mph in 3.1 seconds). The resulting frequency ramp has a 
slope of 18.12 radians. Lor a slope of A<w = 18.12 rad/s , the PLL would take 450.78 seconds to 
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unlock. Even in an accelerating jet, the frequency ramp caused by acceleration is well within the 
limits of the PEL 

Eor normal operating environments, the PEE can be expected to stay locked for a 
real-world frequency ramp. If an environment where the frequency ramp is too large, or applied 
for too long is encountered, a second PEE with an appropriate center frequency could be used. 

2, Analysis of PLL 

The performance measures of the PEE discussed in the previous section give a good idea 
of the capabilities of the fixed-point PEE, but two further items concerning the PEE need to be 
analyzed: the error at the output due to the fixed-point representation, and why an 8-bit system 
failed. The first will be examined first as it will shed some light on the second. 

a. Errors in the Output of the PLL 

To get an idea of the errors expected at the output of the PEE, an analysis of the 
errors that could possibly accumulate at each stage of the PEE is required. A full discussion of 
error analysis and propagation of errors is beyond the scope of this research; however, to 
determine the size of errors at the output of the PEE, some error analysis is required. This will be 
done by looking at the precision at each stage of the PEE and from that precision a worst-case 
error can be determined. By finding the largest possible error in magnitude at each block, the 
error of the sinusoidal output of the NCO is determined. Errors at the output of each block are 
due to two factors. The first is the error at the input of the block. The second is the error due to 
the fixed-point representation and arithmetic operation of the block. Adding these two error 
results in the error at the output of the block. 

The input signal was assumed to be normalized to amplitude of ±1, thus the best 
scaling of the fixed-point number that includes that range is 5 = 2“'"^ (see Chapter 3 on scaling of 
fixed-point numbers). This results in a precision of h.lOE-S . Recall that precision is defined as 
the distance between two successive representable numbers. If a “round to nearest” 
implementation is used, the largest possible error that can occur is half of this precision or 
3.05E-5 . The convention of = accumulated error will be adopted for the remainder of this 
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section. With this notation, when the real world signal is converted to a fixed-point value, 
e, =3.05£-5. 

For the PD component of the PLL, the input signal is multiplied by the output of 
the NCO. The initial value of the NCO is zero, but it is sinusoidal in nature. The assumption will 
be made here that the initial error in the output of the NCO is due only to the precision of the 
fixed-point number representation. The NCO sinusoid output also has amplitude ±1, thus uses 
the same range and precision. For multiplication, the error in the product is the sum of the error 
in the two multiplicands. The error of both multiplicands is the same, thus after the multiplication 
of the reference signal and the NCO, the error due to the input becomes 
2*3.05£’-5 = 6.1035£’-5. Adding the error due to the fixed-point operation yields a value of 
e^=9.\552E-5. 

The next block adds the in-phase and quadrature-phase components in the PD. 
Again, finding the error due to the inputs and adding this value to the error due to the fixed-point 
operation determines the error. For addition, the error due to the inputs is the sum of the errors. 
This results in a value of 1.831 OF'-4 . Adding the resolution error, the output of the summation 
has a possible error of magnitude = 1.8310F-4-I-3.0517F-5 = 2.1362F-4 . 

The next step of the PD was to divide the in-phase and quadrature-phase signals. 
To determine the error at the output of this operation, assume first that the division operation is 
perfect and that the output error is caused only by the error at the input. This being the case, the 
error propagates the same as multiplication. The resulting error is the sum of the errors at the 
input. Therefore the error at the output due to error at the input is 2 * 2.1362F - 4 = 4.2724F - 4. 
Unfortunately the errors of this operation do not stop here. If the denominator in this operation is 
small compared to the numerator, a large number will result. In fixed-point arithmetic, this means 
a larger range of numbers needs to be represented. For the division in the PD, the range of 
numbers was determined to be between ±128 . To get this type of range, a scaling of 5 = is 
needed, which results in a precision of 3.906F-3 . The worst-case error at this stage is the error 
due to the input plus the error due to fixed-point scaling. Adding these two errors together, 
=2.3804F-3. 


49 



The next bloek of the PD is the tan ' (x) funetion. This function is asymptotically 

bounded by ±nH, thus a scaling of 5 = 2“'"^is used. Because this is a non-linear function, errors 
do not work exactly the same. The error will be maximized where the function has a maximum 
derivative. This is at x = 0 . To find the error at the output, the function must be applied to the 
error at x = 0 . More specifically, if the error at the input of this stage is =2.3804£'-3, the 

represented value of the fixed-point system could be 0 or 2.3804£’-3 . If the tan”' function in 
the system was perfectly continuous, the error at the output could be as large as 

=tan”'(0)-tan”'(2.3804£-3) 

= 2.3804£-3. 

Including the error due to the fixed-point representation of tan”' function, the error is 
= 2.3804£'-3-l-3.0517£'-5 = 2.4109£'-3 . As this is the last block of the PD, the total error 
of the phase detector component is e^ pj^ = 2.4109£’ - 3. 

The loop filter component is the biggest contributor to error in the system. This is 
because the loop fdter lead-lag constants are large, thus a larger range is needed, reducing 
precision. The input to the loop filter is the output of the tan”' function, which has a range of 
±71H. This signal is multiplied by the constants ho and hj. For the PLL model, these two values 
were 635.14± 0.015625 and -621.43 ±0.015625, respectively. The output of this multiplication 
needed to include the entire range of ±(.;r/2)(635.14) ~ ±997 . This required a scaling 5 = 2”% 
which has a precision of 0.03125 or output error of 0.015625. Adding the error due to the inputs 
of the multiplication to the error from fixed-point multiplication, the accumulated error becomes 
= 3.366IF'-2 on each multiplication. 

The error gets worse when the addition in Equation (4.5) of the loop filter is 
applied. This addition has three terms in it. The first term in Equation (4.5) has the same range as 
the output of this addition. Because of the feedback of the system, a range of ±5000 is needed. 
This results in a precision of 0.25 or error of 0.125. The error in the three terms are added 
together then added to the precision of the output. The total error of the loop filter then becomes, 

=2*3.3661F-2 + 0.125 + 0.125 = 0.2837. (4.11) 
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The first block in the NCO performs the operation 

OXn) = 0)QTs+Uf{n)Ts. (4.12) 

The first term in this sum is a constant. Evaluated to be 0.62831, it is the center frequency of the 
PEL multiplied by the inverse of the sampling frequency. Eor this PEE model 
=1/Ej =5E-5 . The second term varies with the signal Uj {n); hence a range of numbers is 

required vice the constant first term. The range of u^in) is the range of the output of the loop 
fdter, which was previously determined to be ±5000. Thus the range of the second term in the 
sum is ±5000=±0.25 . Eor this small range, a scaling of S = 2“'^ is used, which increases 
precision. To determine the error at the output of this multiplication, the error due to the precision 
of the output is added to the largest possible error at the input, multiplied by . Thus the 

error at the output becomes, 

+i2-'-’=2.9442£-5. 

Adding the first term in the sum of Equation (4.12), a new range is needed due to the larger value 
of 0.62831. This larger range results in a scaling of 5 = . The error due to this precision is 

3.0517£'-5. Adding the error in the fixed-point representation of 0.62831, the error in the 
second term of 2.9442£'-5 and the error due to the fixed-point operation yields an error of 
e, =9.0476E-5. 

The final possible introduction of errors in the system is due to feedback property 
of the NCO and the continual scaling of the output O^in)io keep it between -7t and.7r. The 
output phase of the NCO equals, 9^{n) = 9j^{n-\) + C 0 QT^+Uf{n-\)T^. If is assumed to 

remain between -7t and 7t , 92{n-\) must be in this range as well. The appropriate scaling for 

this sum is 5 = 2“'^. This results in an error of 6.1035E' - 5 . The last two terms of the sum have 
a total error of 9.0476£'-5. Adding the errors due to fixed-point scaling and errors due to the 
feedback input, the resulting error of ^ 2 (^) is equal to 
=9.04756E-5 + 2*6.1035E-5 = 2.1254E-4. 
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To keep bounded, the algorithm adopted was to subtraet In whenever it got 

larger than n. To get an aceurate error of the system, this subtraetion must be taken into account, 
despite the fact that it occurs infrequently. To subtract In from the phase whenever it grows 
larger than n requires converting In to a fixed-point value. The error in this conversion is 
6.1035£-5. Again adding the error due to the inputs to the error due to the fixed-point 
operation, the total error of the NCO is determined to be, 

3 •3462£-4. (4.13). 

The error determined thus far is the error in the phase of the NCO. The output of 
the NCO is a sinusoid (actually two sinusoids, the in-phase and quadrature-phase components of 
the signal, but the error will be the same for both) with amplitude one. The scaling used for this 
range is 5 = 2“'"^. To determine the output error of this sinusoid, the error at the input must be 
applied to the sinusoid function where its derivative is maximum. For the sine function, this is at 
zero, thus the error at the output due to the input error is sm(3.3462£'-4)-sm(0) = 3.3462£'-4 . 

This value must be added to the error of the sinusoid due to its precision. With a scaling of 2“'"^ a 
fixed-point quantization error of 3.0517£’-5 is found. Adding these two errors together, the 
total error of the system is determined to be 

C,TOr,i=3.3462£-4 + 3.0517£-5 = 3.6513£-4. (4.14) 

The subscript (1) refers to the error after the first iteration of the feedback loop. At 
the beginning of this analysis, the NCO error was multiplied by the reference signal input. 
Because this is a feedback loop and a start point had to be determined, the NCO input error was 
assumed to be zero. Now that a figure for the error in the NCO has been found, this assumption 
requires correction. If the error in the NCO signal is now added and the whole error analysis 
repeated from the PD, the error output after the second iteration is, 

C,ror,2 = 4.5743^-4 (4.15). 

This error analysis shows that as the feedback system continues to operate, the 
error can grow. At this point an idea of how fast errors can propagate and grow in the feedback 
system is determined. By defining. 


C,i C,i-i ’ 
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(4.16) 



the for z = 1 is = 9.2296E -5 . This does not mean that the error will grow at this rate, 
because errors can cancel themselves out. However this is the maximum rate at which they could 
be expected to grow. 

Figure 15 shows the output of the fixed-point NCO subtracted from the output of 
the corresponding floating-point NCO. This difference represents the error in the fixed-point 
implementation. The error can be seen to oscillate, with a maximum value of around 2AE-3 . 
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Figure 15. Error Plot of the Output of the Fixed-point NCO versus Time 


How does this affect the performance of the PEL? Eor the SNR, = 15dB used for 

this simulation, it essentially does not. Using a maximum error of 2.4E' - 3, and the input SNR of 
15 dB, the error reduces this by 0.164 dB. The new . This is not a problem for the PEE 
considered. The PEE locks consistently and as shown in the sections above, within expected 
theoretical limits. The quantization error of a fixed-point system could be a concern if 
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implemented in a low SNR system, but for communieation purposes, with SNRs over 10 dB, the 
reduction will not impact the performance of the PLL. 


b. The 8-bit PLL 

An 8-bit PLL was attempted, but as can be surmised, the 8 bits did not give a large 
enough range or a high enough resolution. When one parameter was fit to the dynamics of the 
PLL, the other was too small, causing critical errors. This last section examines why this was the 
case. 

The only component of the PLL that could not be made to fit an 8-bit 
representation is the loop filter. This is because in this component the largest ranges are 
encountered. Recall from Equation (4.5), the output of the loop fdter is described by 

Uj(n) = Uj^(n-\) + (n) -I- b^u^(n-\) , 

where b^ = 635.1475 and hj =-621.4397. 

The initial cycle of the PLL sets m^(0) = m^( 0) = 0. This means 
Uj{\) = 635.1475m^( 1) . Assuming the worst possible phase offset between the reference and the 
signal sent from the NCO of ^/2, Uj{l) = 7t/2-1.5. From Equation (4.5), 
Uj (1) = 635.1475 *1.5- 982. If a phase offset of -TT / 2 is used, (1) - -982, thus at a minimum 

this range must be included. To achieve this large range with 8-bits requires a scaling of S = 2^. 
The corresponding precision for this scaling is 8. The second iteration comes along and, because 
a correction has been made to the NCO output, the phase error is slightly less, around 1.4. 
Plugging in the formula for Uj , 

(2) = 982+ 933-932 = 939. 

ft is important to notice that the last two terms of this sum add to 1, but with a 
precision of 8, the result is zero, thus m^( 2) = 982 . No correction from the second iteration is 

made. The PEE cannot lock as the corrections to the phase of the NCO frequency cannot be 
realized due to the large range, coupled with the need for high precision. 


54 



Numerous methods were tested to correet this error, but with little sueeess. The 
sum was broken down into numerous sums, so that a higher precision could be obtained, but the 
initial loop made this approach unrealizable. A second more promising approach was to use two 
8-bit numbers in the loop filter, similar to a double floating-point number. Unfortunately, 
Simulink does not make it easy to implement this type of representation and the resulting circuitry 
and eventual hardware implementation complicated the system more than using 16-bit circuitry. 

c. Effects of Changing PEL Parameters 

If not already realized, it should be mentioned that the scaling of the fixed-point 
numbers is a crucial step in developing a fixed-point model. The range of the numbers and 
precision must be considered and evaluated at each step, and when the two conflict, priorities 
must be set to determine a solution. 

It must be realized that changing parameters can result in changing the range or 
precision of certain points in the model. This is especially true when changing the sampling 
frequency of the system. The sampling frequency can have a large range in itself and thus a 
dramatic effect on the PLL operation. For the model considered here, a sampling frequency of 20 
kHz was used. Reducing the sampling frequency to 14 kHz, the model reluctantly locks, but 
jumps are seen at the output, indicating unlock at certain places due to out of range errors. Drop 
the sampling frequency down to 12 kHz and the model is nearly useless as a PLL. The PLL does 
seem to lock, but at the steady state, is unlocked as often as locked. The output signal of the 
system is extremely jumpy, and useless in a communications application. 

Changing the lock range of the system can also have dramatic effects on the 
system. Once the lock range is given, the natural frequency is calculated from the given lock 
range. This value is used to calculate nearly all of the parameters in the PLL. Changing the lock 
frequency can cause an avalanche of changing values, resulting in a different scaling requirement 
for the fixed-point representation of values. 

E, CONVERSION TO HARDWARE 

To complete the analysis of the fixed-point PLL simulation, converting this model to a 
hardware implementation, either on a DSP or FPGA chip needs to be mentioned. Mathworks has 
a product that works in conjunction with Simulink to produce C/C++ code from Simulink models. 
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It provides a debugging feature and eode optimization paekage. This software paekage is ealled 
Real-Time Workshop. At the time of this writing, the software paekage was not available to the 
author, thus testing the eonversion to eode was not possible. The literature suggests that this 
paekage will translate a Simulink model, whether diserete, eontinuous, or fixed-point, into a 
portable, working program. 

Onee this program, in C/C++ or Ada, has been developed, eonverting it to a VHDL 
language is simply a matter of syntax (HDL stands for Hardware Deseription Language). This 
type of language was adapted as the need for a language that could be used by a program- 
controlled machine for generation of final hardware was recognized. Several commercial HDL 
packages are available today. 


In this chapter, the procedure for designing the fixed-point PLL model was described. 
This was done by building a floating-point model and converting it to a fixed-point equivalent. 
Analysis was done by comparing the fixed-point simulation results to those obtained from the 
floating-point. Finally an error analysis was done to determine how large errors could be 
expected to grow due to the fixed-point arithmetic quantization error. The final chapter deals with 
capabilities, limitations, and future research for a FPGA based PLL. 
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V. CONCLUSIONS 


The objectives of this research were to build a fixed-point PLL and analyze the 
performance of the model. The was accomplished by building a floating-point PLL model and 
converting it to fixed-point arithmetic and comparing the results of the two simulations. This 
Chapter discusses the capabilities and limitations of the fixed-point model, and concludes with 
recommendations for further research in this area. 

A. CAPABILITIES AND LIMITATIONS 

The fixed-point model developed was able to lock onto a given reference input signal with 
comparable performance to the floating-point or analog equivalent. Lock-ranges and lock-times 
in the model were equivalent. The model proved to be able to track a drifting frequency, lock 
onto a frequency step larger than the theoretical value, and a phase step with accuracy and within 
theoretical limits. As a fixed-point implementation, it is a viable option for use on a FPGA or 
DSP chip, fulfilling the requirement for a PLL in a software radio package. The reason this is the 
case is because the fixed-point operations are much faster than the IEEE floating point standard, 
allowing real-time radio in a downloadable software-like format. The second major advantage of 
the fixed-point implementation is the memory storage required. Because the range of the 
numbers is known, 16-bits can be used and scaled appropriately to represent each number as 
opposed to 32-bits or 64-bits of floating point. This reduces the amount of memory required by 
50% or 66%, a valuable commodity on a small DSP chip. 

The cost of fast operation and memory savings is a trade-off in precision and range. Using 
only 16-bits means the range of numbers or the precision much be sacrificed. The limitation of 
the fixed-point model is the parameter specific operation of a given model. The point behind 
using 16-bits and a fixed-point operation is that the range of numbers is known, and the window 
of numbers that 16-bits can represent can slide to accommodate that range. The other side of this 
is that, if parameters are changed and that window requiring representation changes, the fixed- 
point model could fail. To have a completely portable model, the scaling of the fixed-point 
numbers needs to be variable. By allowing the scaling to be variable, any change in parameters 
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can be accounted for by changing the scaling, hence sliding the window of fixed-point 
representation to match the changing range. 

The 16-bit fixed-point model introduced error into the PLL that otherwise would not be 
there. This is seen in the quantization error in representing a number. To speed up arithmetic 
operation, fixed-point numbers are represented in a way that makes binary arithmetic - the 
arithmetic a processor uses - much faster. This representation is not optimal for representing 
real-world values. As was shown in Chapter 4, this error resulted in a reduction of the SNR of 
around 0.16 dB. For radio communication purposes, this is not an obstacle, but for applications 
using a small SNR, limitation on fixed-point representation needs to be considered. 

As mentioned numerous times, fixed-point implementation has a smaller range of 
representable numbers, given a certain precision. For the loop filter used, a certain precision is 
required to make sure small additions have an effect on the output. This results in the precision 
being the priority and, hence, the range suffers. Two limitations are the cause of this limited 
range. 

First, an 8-bit PLL was not realizable. By using 8-bits, the range was decreased to the 
point that numbers were over 50% outside the representable range. This produced critical failure 
errors in the feedback loop, preventing the PLL from locking. Several solutions were attempted 
and could hold promise, but a 16-bit fixed-point implementation was preferred. 

The second limitation of the reduced range is that the pull-in range of the system was 
reduced. While a second-order analog or floating-point PLL with an active P-I loop filter should 
have unlimited pull-in range, the fixed-point counterpart did not. This was due to the range 
restrictions imposed by the fixed-point representation. Because only a limited range could be 
represented, when a frequency outside that range was encountered, the value could not be 
approximated in the fixed-point scaling. This inevitably prevented the system from locking. 


B, RECOMMENDATIONS FOR FURTHER RESEARCH 

The main area for further research is to actually build a phase lock loop on an FPGA or 
DSP chip. This step will have two steps. First, the PLL model must be converted to useable 
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code. The Mathworks Real Time Workshop is ideal for this purpose. Researeh into this program 
and the how eode is produeed from a Simulink model needs to be done. This researeh would 
require understanding how different parameters in Real Time Workshop affeet the generated 
eode, and how to optimize the final program. 

The seeond step would involve translation to an HDL . While there are several 
eommereially available HDL’s, ineluding Interaetive Design Language (IDL) from IBM, 
Instruetion Set Proeessor Speeiiieation (ISPS), Test Generation and Simulation (TEGAS), Texas 
Instruments HDL (TI-HDL), and ZEUS, ereated at General Eleetrie Corporation, the most 
promising language was developed in 1983. Contraeted by the DoD, VHDL 2.0 was released in 
late 1983. In 1987, VHDL beeame the IEEE standard for HDL. It is suffieiently rich for 
designing digital systems and supports a hierarehieal description of hardware. This means it ean 
be used to deseribe programs from systems to gates, or even the switeh level [7]. To continue this 
researeh, an understanding of an HDL and a conversion from the present Simulink model needs to 
oeeur. This proeess is not as easy as a simple eonversion from C/C++. Erom the C/C++ eode, a 
Preliminary Design Deseription (PDD) is typieally generated prior to VHDL eoding. Onee an 
aeeeptable PDD is eomplete, VHDL eoding ean eommenee [3]. 

Analysis on what type of LPGA or DSP is needed to satisfy the requirements of a software 
radio needs to be done. This would inelude figures of merit sueh as speed and memory. If the 
radio requires GHz operation, but an EPGA ean only handle MHz, a real-time system is not 
feasible. This is not the ease, as “IE Sampling” ean bring the operational frequeney down, but 
researeh on the trade-off between high speed and affordability is required. Researeh into how 
much memory would be required to not only put the PEL on the ehip, but the other eomponent 
parts of a radio sueh as a demodulator, deeoder, or delay loek loop (DLL). 

C. FINAL COMMENTS 

The final goal of this and related researeh is to build a software radio on a DSP or 
preferable an EPGA ehip. The idealized software radio would allow the user the flexibility to 
download software paekages related to their mission at that time. This platform would be a 
general-purpose radio. Instead of using several radios for different situations, modulation 
formats, or radio protoeol, a single unit, eapable of eommunieating with all formats by a 
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downloadable software package would be used. By doing all the demodulation, and signal 
processing in software, this goal can be achieved. A typical radio signal would be digitized at the 
antenna, or more likely, an analog frequency translator would be used to shift the RF signal to an 
intermediate frequency to reduce the sample rate requirement of the A/D converter. With the 
signal digitized, all radio processing function are done in software. An integral part of this radio 
function is the PLL. To realize these functions on an FPGA chip, fixed-point arithmetic is 
required. For reasons already specified such as memory space and speed, DSP and FPGA chips 
use fixed-point representation rather than floating point. Thus a fixed-point PLL model was 
developed. 

The fixed-point PLL model simulated in this research performed comparable to an analog 
or floating-point equivalent. The fact that it was implemented using fixed-point arithmetic makes 
it ideal for conversion to an FPGA or DSP chip. 

The limitations introduced by the fixed-point arithmetic must be considered in designing 
the final product. Understanding the nature of fixed-point scaling and how parameters can affect 
the ranges of an application are essential. For the model in this research, a 16-bit fixed-point 
representation was assumed. An 8-bit model would have been ideal, but the limitation on range 
and precision of only 8-bits caused a deferment to 16-bits. Fortunately, DSP and FPGA chips 
come in 16-bit and even as high as 32-bit on the commercial market, making the 16-bit 
implementation a practical solution. 
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APPENDIX A 


This Appendix gives the eode for the MATLAB floating-point PLL and the fixed-point 
Simulink model used for simulation. 

Floating Point MATLAB code 

For ease of referenee, the variables used in the eode will be matehed to the notation used 
in the text. 

• sigi = in-phase eomponent of w,. 

• sigQ = quadrature-phase component of Wj. 

• zeta = ^ . 

• Wl= co^, Wn= co^. 

• taul = fj, tau2 = ^ 2 . 

• phi = 02 • 

• theta = 9^. 

• FI = output frequency of U 2 . 

The following is the code generated for PLL model. 
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ADPLL2 (All-Digital Phase Lock Loop) 

Implements a phase lock loop assuming a complex in-phase and 
quadrature-phase input signal. Adapted from Reference [1]. 

Michael Johannes 
12 September, 2002 
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% Initialize signal 

parameters 

Fc = 

2100; 

g, 

o 

Input Reference signal frequency 

UlO = 

1; 

g, 

o 

Input Reference signal amplitude 

U2 0 = 

1; 

g, 

o 

Output NCO signal amplitude 

FO = 

2000; 

g, 

o 

Center Frequency of NCO 

Fs = 

40000 

g, 

o 

Sampling freqency 

Ts = 

1/Fs; 

g, 

o 

Sampling duration 

% Construction of 

Reference signal sigi and sigQ 

SNRdb 

= 15; 

g, 

o 

Input Signal to Noise ratio in dB 

BW = 

1000; 

g, 

o 

Prefilter bandwidth 

pod = 

rand*3 60; 

g, 

o 

Randomize phase offset in degrees 

por = 

pod*pi/180; 

g, 

o 

Convert random phase offset to rad/sec 


t = 0:Ts:50/Fc; % Vector of time samples 

% Calculate noise power from SNR and prefilter bandwidth 

sigma = sqrt(.5/(10^(SNRdb/10))/BW) 
sigma = sqrt(.5/10^(SNRdb/10)/BW); 

% Signal construction with noise added 

sigi = U10*cos(2*pi*Fc*t+pi); 
sigQ = U10*sin(2*pi*Fc*t+pi); 
sigI = sigl + sigma*randn(size (t)); 
sigQ = sigQ+sigma*randn(size(t)); 

% Calculation of PLL parameters 


zeta = .6; 

g, 

o 

Damping factor 


W1 = 200*pi; 

g, 

O 

Lock range defined according 

to expected 


g, 

o 

input characteristics 


Wn = Wl/2/zeta 

g, 

O 

Natural frequency calculated 

from lock range 

Kd = 1; 

g, 

o 

Phase detector gain 


Ko = 1; 

g, 

o 

Loop filter gain 



% Loop filter constants and lead-lag gain constants calculated 
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taul 

= Ko 

*Kd/(Wn^2 

tau2 

= 2* 

zeta/Wn; 

al = 

-1; 


bO = 

(Ts/ 

(2*taul)) 

bl = 

(Ts/ 

(2*taul)) 

% In 

itialization 

phi = 

0; 


theta 

= 0 

r 

Uf=0; 



theta 

I = 

0; 

thetaQ = 

0; 

FI = 

2*pi 

*F0 

% General 

case pha 

for n 

= 2 

:length(s 


% Lead constant 
% Lag constant 


Output phase of the NCO 

Phase error of the two signals 

Output signal of the loop filter 

Error of in-phase signal 

Error in quadrature-phase signal 

Output frequency of NCO 


gain 

gain 


phi (n) = phi (n-1)+F1 (n-1)*Ts; % Phase error from current cycle 

% added to the phase of the 

% output signal 

% Subtract 2*pi whenever the phase is larger than pi to keep the 
% phase bounded 

if phi (n) > pi 

phi(n) = phi (n) - 2*pi; 

end 


% Implementation of phase detector 

thetal (n) = sigi(n)*cos(phi(n) )+sigQ(n)*sin(phi(n) ) ; 
thetaQ (n) = sigI(n)*sin(phi(n) )-sigQ(n)*cos(phi(n) ) ; 
theta(n) = atan(-thetaQ(n)/thetal (n) ); 

% Loop filter equation. 

% Implements the transfer function l+s*taul 

o, 

o 

% s*tau2 


Uf(n) = -al*Uf(n-1)+bO*theta(n)+bl*theta (n-1); 

% Implements the equation for the NCO Wo+Uf(n) 
Fl(n) = 2*pi*F0+Uf (n); 

end % for loop 
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Fixed-point Simulink Model 

The Simulink model was broken down into its eomponent parts; Phase Detector, Loop 
Filter, and Numerically Controlled Oscillator. Figure A1 shows all components in the PLL 
system. Figures A2-A4 show the Simulink block structure for each component. The scaling of 
the fixed-point representation for each block is shown next to that block at the lowest level. 



Figure A1. Phase-Lock Loop Simuli nk Model 

The PD takes three inputs. The inputs Ini and In2 are the in-phase and quadrature-phase 
reference signals, depicted by the [t', sigl'] and [t', sigQ'] , respectively. These signals are 

generated in the MATLAB workspace, and imported into Simulink using the From Workspace 
blocks. The t input of that block simply tells the model it is a sampled signal and the duration 
between samples. The input In2 is the phase error feedback from the NCO. 

The FixPtGateway Out blocks take the fixed-point value and convert it back to a double 
floating-point representation so that it can be graphed versus time in the Scope block. 

Figure A2 is the phase detector block structure embedded in Figure Al. The three input 
lines are seen at the left of the figure. Inputs Ini and In3 are the sampled signal and immediately 
go into a FixPtGateway In block that converts the floating-point values to a fixed-point 
representation. The final input, In2, is the total phase of the NCO. It serves as input to a block 
labeled sin/cos. This block takes the phase and by use of a look-up table, outputs the sine and 
cosine of the input phase. The output, Outl, of the PD is the phase error of the system plus higher 
frequency terms. 
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Figure A2. Phase Deteetor Simulink Model 


The next block from Figure A1 is the loop filter. This block is shown in detail in Figure 
A3. The constants 635.146 and -621.44 are the lead-lag constants andbj. The 1/z blocks 

signify a unit time delay. This block saves the output from the last sample time and outputs it at 
the current time sample. The output of this component is the phase error of the system. 
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Figure A3. Simulink Active P-I Loop Filter 

The final component block of Figure A1 is the NCO. This subsystem is shown in Figure 
A4. It takes as input the error from the loop filter and calculates the total phase error of the 
system. This phase error is fed back to the phase detector, and the process repeats until the error 
is forced to zero. The complicated circuitry at the end of the NCO is a comparator and switch 
combination that compares the phase output io n. If the output is less than n , the switch outputs 
its top input, which is just the signal. If the phase is greater than Tt , the switch outputs the bottom 
input, which is the signal minus iTt. This circuitry keeps the phase from growing unbounded or 
out of range of the fixed-point model. 
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Figure A4. Simulink Numerically Controlled Oscillator 
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